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ABSTRACT 


This thesis is a tutorial on digital signal sodulation 
techniques used in satellite communications and includes 
computer simulations of those digital signal sodulation 
technigues introduced. The purpose of the thesis is to 
introduce digital signal modulation technigues ard through 
the use of computer simulation, generate statistics which 
represent the characteristics of the FFT for the respective 
Signal type. Further, an analysis of tne statistics of the 
FFT's was conducted to determine if there is any relation- 
Ship between the components of the FFT of the different 
signals. The statistic used to investigate this possible 
relationship was the F-distribution. The computer simula- 
tion was written and conducted in the FORTRAN programming 
language. A copy of the program, results of the simulations 
and the statistical analysis conducted are included ‘in the 
appendices. , 
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A. BACKGROUND 


Since the introduction of the sampling theorem and tae 
matched fiiter, digital communications techniques have 
developed into a highly proficient discipline. The marriage 
of this discipline with the rapidly expanding space frogran 
has rcesulted in ccmmunication satellites employing a 
multitude of digital signal modulation tecnnigques. A 
modulation technique is a method of transmitting the 
information contained “in a message by varying cr modulating 
the characteristics of a carrier waveforn. These metnods 
offer a range of advantages and disadvantages depending on 
the specific characteristics of the modulation technique 
employed. ‘The applications of the various signal moduiation 
tecanigues likewise vary. In order to understand this 
exciting new field it is necessary to look at some of the 
aspects of satellite communications systems in general. 
Then an investigation will be made into the general 
attributes of digital communications and their relationsnip 
to the satellite system 4 , 

The ability to understand these digital signal 
modulation techniques is the first step in being able to 
intercept, identify and demodulate unknown digital signals. 
These digital signals, transmitted from unknown sources, are 
believed to display frequency characteristics pecuiiar to 
the mcdulation technique employed in the encoding process. 
The digital computer offers unique opportunities in 
Simulating these aodulation techniques, in siynal frocessiag 


and in decoding of an intercepted signal. 
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B. SPECIFIC GOALS 


It is the specific goal of this thesis to investigate 
the open literature on the topic of digital signal 
modulation technigues in satellite communications. This 
includes a basic understanding of the communications 
satellite system as well as the specific techniques employed 
in signal nodulation. Additionally, once a basic 
understanding of these digital signal modulation techniques 
is achieved, computer code in the FORTRAN programming 
language will be developed which simulates these modulared 
Signals. The statistics of the time-varying Fast Fourier 
Transforms (FFT) of these Simulated signals will tLe 
investigated. The purpose of this analysis will be to iead 
to fcllow-on research in the area of signal analysis of 
intercepted digital signals whereby they can be classified 
by their FFT as employing a specific digital signal 
modulation technigue. 


C. SCOPE OF THE PROJECT 


Althouga an indepth analysis of all the digital signal 
modulation techniques which will be introduced involves 
considerable higher level mathematics, it is not-within the 
scope of this project to delve heavily into the mathematics. 
It would be advantageous, however, for the reader to have 
had integral and differential calculus and an introduction 
to statistics. Also a course in electrical engineering may 
prove helpful but is not essential since in actuality it is 
the electrical engineering aspects of satellite 
communications that this tutorial is attempting to present. 
The computer program@ing will be accomplished utilizing the 
techniques of software engineering and top down modular 
design. existing blocks of code will be used as modules 
whenever appropriate and available. Again it is the 
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ultimate purpose of this project to examine a variety of 
digital signal modulation techniques and develop computer 
code that simulates common signals used in satellite 
communications that have been produced by one of the many 
digital signal modulation techniques to be investigated. 
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Il. JBDERSTANDING SATELLITE CONNONTCATIONS 


A. HISTORY AND APPLICATIONS 


Forty years ago, in 1945, Arthur Clarke first envisioned 
the use of space stations placed in geosynchronous orbit for 
communicating to different points on the earth [Ref. 1]. 
Nine years later in 1954, J.R. Pierce of Bell Lakoratories 
performed a system analysis on such a communications systea 
(Ref. 2). In 1957, the launch of Sputnik demonstrated the 
feasibility of using a satellite for just such an 
application. However, by 1961 the only satellite 
communications technologies which had been demonstrated were 
the Ccurier 1B satellite, a short-life active retransmission 
teletype communications satellite ina 1000 km orbit, aud 
the Echo I balloon in a 1600 km orbit which demonstrated 


- passive reflection of powerful microwave signals from one 


earth station to another. Active microwave communications 
demonstrated in Projects Telstar and Relay were years away 
(Ref. 3]. The first geostationary orbit was achieved by 
Project Syncom in 1963 [Ref. 4}. 

In 1964, communication organizations from several 
countries joined together to form the international 
organization of INTELSAT (International Telecommunications 
Satellite Organization). . INTELSAT's purpose was to develop 
a satellite network which would provide truly glotal 
communications capabilities. This resulted in the launch in 
1965 of the world's first commercial communicatioas 
satellite, INTELSAT I, also known as “Early Bird". With 
“Early Bird", telecommunications utilizing satellite relay 
were established between the United States and Europe. 
(Ref. 3] . 
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The reliability of these satellite systems has improved 
dramatically since INTELSAT I in 1965. Reliability of 
individual links in the system approach 99.99 percent or 
higher. Total system reliability exceeding 99.9 percent is 
common [Ref. 5], naking Satellite communications more 
reliable than gany other modes of communications. In this 
sense, reliability is a measure of the probability that no 
failure will occur in a respective channel or in the systen 
during the design life of the satellite. 

In crder for the operability, capability and reliability 
of these systems to have developed at this pace, it was 
necessary for the technologies associated with them to 
develcp as rapidly or even more rapidly than the systems 
themselves. Engineering sciences and specifically the 
fields of aerospace and electrical engineering historically 
have required between 7 and 10 years to take a concept from 
operational requirement to full scale operation. This was 
not the case with the concept of communications satellites 
which das seen four generations ot INTELSAT sateiiites 
within a decade's time. [Ref. 3} 

The "Early sBird" satellite weighed approximately 32 
kilograms and possessed limited power and bandwidth capacity 
enabling it to carry only 240 two-way telephone 
conversations. Today's communications satellites represent 
order-of-magnitude improvements in many important operating 
parameters such as fower and bandwidth. Additionally, 
increased effective radiated power from the techniques of 
stabilized earth pointing antennas have greatly increased 
the capacity of later generation communications satellites. 
INTELSAT V, the current generation of communications 
satellites, is a three-axis stabilized platform using not 
Only the 6/4 GHz frequency band (6 GHz receive/4 GHz 
transmit) as in earlier generation satellites but also a 500 
MHz kandwidth available in the 14/11 GHz band. The 
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separation in the receive and transmit frequencies is 
necessary to prevent interference during simultaneous 
operation ot the receiver and transmitter. Ancther factor 
increasing the capacity of present generation satellites is 


the increase in prigary power available in the satellite. 
(Bef. 3) 
The technologies which have contributed to the evolution 


of coggunications satellites come from two primary sources, 
namely technology from the space program of the 1960's. 
supported by NASA and DoD and communications technology due 
laryely to commercial and private sector contributions. 
(Ref. 3) 

Electronic devices and ‘gmponents have contributed 
significantly to the rapid development of satellite systems. 
These electronic devices range from something which is now 
considered bkhasic, i.e., the transistor, to devices such as 
Traveling Wave Tube aaplifiers, lightweight antennas and 
antenna feed systems. There have been najor improvements in 
satellite power sources including more efficient solar ceils 
and high storage capacity/lightweight batteries. 
Additionally, two Significant aevelopments in electrical 
engineering have made modern day digital signal processiny a 
reality. They are the sampling theorem and the natched 
filter. (Ref. 3] 

Communications satellites of the future are likely to 
utilize onboard signal processing. Sigral processing 
functions such as signal reshaping, Switching and/or 
multiplexing and compression could soon take place on the 
satellite due to the reducticn in size and weight of the 
necessary hardware. Operation at over 100 megabits per 
second are envisioned. Onboard signal processing will 
greatly reduce the exrense and complexity of earth stations 
thereby making services of a satellite available to a wider 
range of small and gecgraphically disperse users. [{Ref. 3] 
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Satellite comsunications . hold great promise to provide 
service to a multitude of users over a wide area. 
Applications lie not solely in retransmission but also in 
data collection from that same large geographic area. 
Military users have definite applications in these areas 
when warning, intelligence and surveillance systems provide 
digital inputs to a central source. Although most of the 
present applications for satellite communications are still 
provided by the government, commercial applications nave 
seen tremendous increases in the last several years. 
(Ref. 3] 


Be. ORBITS AND LIMITATIONS 


Satellite orbits are generaily categorized as either 
equatcrial (0 degrees inclination), polar (90 degrees 
inclination) or inclined at some angle other than 0 or 90 
degrees relative to the spin axis of the earth as 
illustrated in Figure 2.1 (Ref. 6]. Each satellite orbit 
has a characteristic velocity which is dependent on the 
height of the orbit and the orbit's eccentricity. 
Eccentricity is a measure of the degree to which the orbit 
approximates a circle, A circle has eccentricity equal to 
zero since the focus is located at the center. See equation 
2.1 and Figure 2.2. 

A communications satellite in elliptical orbit about tne 
earth obeys Kepler's second law of planetary motion. That 
is, a satellite's constant angular momentum about the earth 
means that its areal velocity also remains corstart. See 
Figure 2.3. Of particular interest is the satellite 
velocity at apogee (Va) and perigee (Vp) given by equations 
2.2 and 2.3. 

Note that for a circle, e = 0 and Ra = Rp. Therefore a 
satellite in circular orbit atout the earth has an orbital 


velocity as given in equation 2.4. 
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Figure 2.1 Satellite Orbits 


Most communications satellites are placed in a circular 
equatorial orbit where the desire is to stabilize them over 
a fixed point on the surface of the earth called the 
subsatellite point. This type of orbit is called 
geostationary or stationary. Any satellite which has an 
orbital period equal to the period of rotation of the earth 
is called synchronous or geosynchronous. These terms are 
used almost interchangeably in most literature. AS 
mentioned for a geosynchronous conmunications satellite, the 
orbital period of the satellite, T, must be equal to the 
period of rotation of the earth. Kepler's third law of 
planetary motion can be rewritten to yield equation 2.5. 

The period of revolution of the earth for a Sidereal day 
is 23 hr 56 min 4 sec. For that given period there is only 
one satellite altitude as expressed by Kepler's third law. 
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Figure 2.2 Eccentricity 


eccentricity (e) = c/a (eqn 2.1) 


distance from focus 


semi-major axis 


b = semi-minor axis 
Ra 
Rp 


radius cf apogee 


radius cf perigee 


By vtearranging terms that altitude is given in equation 2.0. 
Since the orbit is circular, a= Ra = Kp and the héigat of 
the orbit above the surface of the earth is o = a-=—- Re; or 
35,804 kn. 

One disadvantage of a geosynchronous communicationos 
satellite is the lack of global coverage. These satellites 
provide excellent coverage of the subsatellite point anéd 
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no Figure 2.3 Kepler's Second Law of Planetary Motion 
. | ; 
ey Va = (k/Ra (i-e) )**.5 (eyn 2.2) 
he 
te Vp = (k/Rp(1te))**.5 | (eqn 2.3) 
ae k = GMe = 3.98866 x 10**-11 m3/s2 
ge G = universal gravitational constant 
Bae G = 6.67 x 10**-11 N m2/kg2 
' N = newtons = m kg/s2 
Aa Me = mass of the earth = 5.98 x 10**24 kg 
rad | 
$s. Vo = (k/Bc) **.5 (eqn 2.4) 
ay 
am Ro = Re + b 
Re = radius of the earth = 6.37 x 10**6 
os h = satellite orbital altitude 
a 
pe T = 2 pi a¥#1.5 / k¥*#.5 (eqn 2.5) 
ue = period of rotation 
van = semi-major axis 


AT / 2 pi) **2/3 (k) ** 1/3 (eqn 2.6) 


o 
" 


42,173 ko 
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laterally to a latitude of about +80° {Ref. 6]j. This is 
generally no problem for commercial applications since the 
polar regions do not require a significant degree of access. 
The military, on the other hand, does have an interest in 
communications in the polar region and therefore has several 
Communications satellites that have orbits inclined at 
various angles. This type of orbit generally has 
disadvantages of lack of continuous coverage and a much more 
complicated systen of ground tracking and receiving 


staticns. 


C. FREQUENCY BAND CONSIDERATIONS 


Although there appears to be an infinite number of 
frequencies available for communications, restrictions do 
exist as to. those which are practicable. Limitations on 
available frequency bands for satellite communications are 
due tc the need to select segmerts of the electromagnetic 
spectrum which reduce noise and interference and are nost 
favorable in terms of power efficiency and propagation 
distortions. Trade offs must be made to arrive at the 
optimum frequency for a particular application since single 
frequencies seldom cffer the best performance for ali 
variables. The problems arise when consideration is given 
to the number of users reguicring the same frequency bands 
inciuding terrestrial communications networks. Since the 
problem of iaterference is glotal, a worldwide organization 
has been established to assign frequency bands for various 
applications. This organization is called WARC, World 
Administrative kadio Conterence (Ref. 6]. Illustrations of 
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the pertions of the electromagnetic spectrum in guestion and 
current allocations of satellite frequency bands are shown 
in Figures 2.4 and 2.5 (Ref. 6]. 
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Figure 2.4 Electromagnetic Spectrua 


In general, as long as the hardware and technology will 
support it, higher frequencies are more desirable and nore 
in demand because they offer higher tneoretical capacity. 
This is due to the fact that only a percentage of the 
carrier frequency is capable or actually transmitting a 
signal of a given bandwidth. Higher frequencies wouid also 
experience less interference with existing land and 
satellite systems. A further discussion of the bandwidths 
available as a function of freguency will be included at a 
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| = downlink, - = crosslink). 


Nae ce a Me ae So ee a a ne A a A SO ce an pe OT: SS ee SD een 


Figure 2.5 Satellite Frequencies 


later point under the topic of actual digital signal 
techniques. 

As an additional note, in order to assist in the ability 
to handle more signals in an obviously linited 
electromagnetic spectrun, several multiple access schetnes 
have teen developed. The mest common techniques involve 
rrequency-division multiple access (FDMA) where the 
allocated satellite frequency band is divided among the 
users into specific uplink and downlink frequencies. Next 
is time-division mnultifle access (TDMA) where the satellite 
frequency band is shared by all users by carefully dividing 


23 


ste 
aT 
a) 


a) 
+ 


the time any one user has access. Finally. there is 
code~-divisiona multiple access (CDMA) which involves 
modulating a specific address code which has been 
superimposed on the signal directly onto the carrier. In 
this manner all users share the satellite frequency band and 
only those receiving stations that can demodulate the 
address code can receive the specific signal. [Ref. 6] 


D. WHY DIGITAL SODU LATION 


In general, digital signal modulation techniques are 
superior to analog signal modulation technigues used in 
satellite communications for the following reasons: 

1. Compatibility with digital computers 

2. Economic advantages 

3. High degree of flexibility 

4%. Less susceptible to interference 

5. Quality of signal independent of transmission 
istance and network makeup 


1. Compatibility with Digital Computers 


Clearly one of the most distinctive advantages of 
digital techniques over analog techniques involves computer 
applications, Properly formatted digital signals can be 
used to represent any analog signal (more on this under the 
discussion of the sampling theorem). Once in the digital 
format these signals can be easily manipulated within the 
digital computer. Arithmetic operations can te applied as 
well as logical operations, The signal can te _ stored 
without alteration or delayed and therefore can be used to 
"simulate" real physical situations. 

The hardware associated with digital circuits is 
free from drift or aging and does not require caliktration. 
Additionally digital circuitry is compatible with present 
day integrated circuit technology allowing a standardized 
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building block construction approach. The operating 
characteristics of these systems employing a digital 
computer can be changed by altering the software rather than 
the hardware as is the case fcr analog systems. Finally, 
the use of digital computer technology allows tine 
aultiplexing. 


2. Flexibility and Economy 


The flexibility and economy of digital satellite 
communications comes from the fact that more and nore 
processing can be done onboard. This allows the uplink and 
downlink to be compietely separated. This regenerative 
Mature makes it possible for low error rates and aigh 
reliability through the use of digital technigues not 
available to analog systems. Because digital signal 
processing or multiplexing is less costly than for analog 
signals, Simpler and cheaper interfaces between earth 
staticrs and terrestrial communications networks are 
possitle. Additionally, there are reduced production costs 
and increased capacity associated with digital circuits. 
(Ref. 7] 


3. Quality and Interference 


The capability of digital systems to regenerate tue 
Signal and allow for multiple switching and Siynai 
processing without degradatvion in signal quality makes the 
digital Signal basically independent of transmission 
distance. Multiple hops from satellite to earth staticn or 
satellite to satellite are possiple without accumulation of 
the ncise characteristic of analog systeas. Additicrally, 
digital systems are capable of operating at a Signal to 
noise ratio of 20 dB to 30 dB as compared to analog systems 
reguiring a much more powerful] signal. [Ref. 8] 
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TIT. [NTRODUCTION FO THE FOURIER TRANSFORM, SANPLING 
THEO BEM, 


AUTOCORBELATION FUNCTION AND THE MATCHED FILTER 


Essential to the understanding of digital signal 
modulation techniques are a few basic tools of the 
electrical engineer and the mathematician. These tools will 
be developed and elaborated on to the extent necessary to 
understand their applicability to the subject of digital 
communications. The description is not meant to be a 
detailed investigation of the respective topics. where 
relevant, the application of the concept being descrited 


wili ke mentioned. 


A. FCURIER TRANSFOR SB 


‘In mathematical terms, voltages can be expressed as 
functions of time or of frequency. It is more commcn to see 
voltages represented as a function of time as in equation 
33:12 


v(t) = A cos(wt) ; (eqn 3..1) 
A = amplitude 

w = angular frequency = 2 pi f 

£ = frequency = 1/T 

T = period 


It iS important tc note that both the time and frequency 
functions are representations of voltage and as such may be 
used interchangeably. The Fourier representation (v(t) } = 
V(f) ais avoltage descriptor in the frequency domain (a 
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function of frequency) while v(t) is a voltage descriptor in 
the time domain (a function of time). They are different 
but equivalent and either voltage descriptor may be used, 
depending on the concept being explored, to best represent 
the voltage within context. 

The Fourier transform is of principle interest rather 
than the Fourier series since the latter is applicable only 
to periodic voltages. The Fourier transform is applicable 
to voltage pulses, random voltages and other non-periodic 
voltages. 

DEFINITION: 

foe) 

[vty e**-jut dt (eqn 3.2) 
- 0 


Fev(t)] = Vw) 


V(t) <qcennnena-=> VW) 


time domain frequency domain 


Rememktering Euler's formula 
e**-jwt = cos(wt) - j sin (wt) (eqn 3.3) 
the Fceurier transform becomes 


[o ¢) 
Al v(ty] = vw) = f v(t) [cos (wt) - j sin(wt)] (eqn 3.4) 


rove) —0O 
Vv (Ww) = [ v(t) cos (wt) = jf v(t)sin (wt) (eqn 3.5) 
—co — 00 


Since it will be seen that digital communications deals 
priaarily with pulses of finite duration (expressed as 
period, T) , it is worthwhile to examine the Fourier 
transform or a pulse cf amplitude A and duration T. 

A; -T/2 < t < Ts2 
let v(t) = 


0; elsewhere 
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Figure 3.1 Square Pulse 


Figure 3.1 is a representation of v(t) or the voltage 
expressed in the time domain. The position of v(t) on the 
t-axis was chosen for convenience of integration but could 
have been situated anywhere on the time line. 


co 


Fev(ty} = vw) = [v(tyer*-jut at (eqn 3.6) 
y ze ‘Ve : ; 
Ve) = Jere jue dt + [a e**—jwt dt + (eqn 3.7) 
7 Osea e 4ut dt hs 
a Wr 
V(w) = fa e**-jut dt (eyn 3.8) 
W 


The integration above may te attacked either head on or 
by sukstituting cos(wt) - j sin(wt) for e**-jwt. The direct 
approach is illustrated due ta tne relative simplicity of 
the integrand. 


YW 
V(w) =fa e+*- jvt dt (eqn 3.9) 
“%, 
T/2 
V(w) = -A/jwle*e-jut] | (eqn 3.10) 
-T/2 
28 
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V(w) -A/jwle**-juT/2 ~ e** jvT/2] (eqn 3.11) 


V(w) A/jule** jwls2 - e€ *-5uwT/2j (eqn 3.12) 


By sukstituting 2pi f = w, the following result is obtained: 


V(£) = A/j2 pi fl[e**j2 pi f£T72 - (eqn 3.13) 
e**-j2 pi f£T/2] 


Vif) = A/j2 pi fl e**j pi £1 - e**-j pi fT] 


Using Euler's formula, i-e., sin 9 = (e**j0-e**-j8) /2j5: 


V(f) = 2jA/j2 pi f[ (e**j pi fT - (eqn 3. 14) 
e**=-5 pi £T)/2j] 


V(f) = A/pi f£{sin( fi FT) } 


Knowing that sin x / x = sinc x 


V(£) 


AT/pi £T[sin(pi fT) ] (eqn 3.15) 


V (£) AT sinc({pi fT) 


The sinc function is common in digital electronics and 
plots as the product of sSin({pi £T) and 1/(pi fT) as int 
Figure 3.2. In Figure 3.2, V(f£) in the frequency domain is 
equivalent to v(t) in the time domain. Note that as the 
pulse, T, gets longer, 1/T gets smaller or the first zero 
crossing of the sinc function occurs ata lower and lower 


freyuency. 
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Figure 3.2 Plot of Sinc Function 


1. Amplitude and Phase Spectrun 


Recall that V(w) can be represented as in equation 
3.16. The cosine term is the real part while the sine tern 
is the imaginary part. By referencing Figure 3.3, a brief 
review of the complex plane .is accomplished and its 
relationship to the amplitude and phase spectrum of a given 
voltage is represented. See eguations 3.17 through 3.20. 


ve) re 
V(w) = [vtrcos (wt) - ifv (t) sin (wt) (eqn 3. 16) 
— oo real ~ COimaginary 


{V(f) 1 is called the amplitude spectrum of the given 
voltage. The amplitude spectrum can also be calculated by 
using the complex conjugate or the Fourier transform ani is 
always positive as shown in equation 3.21. The phase 
spectrun, 6, of the furction in question is represented by 
the arctanf imaginary/real} and is illustrated in Figure 3.5. 
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Figure 3.3 Phasor Diagran 


c= {a2 + b2] (eqn 3.17) 
a + jb = [a2 + b2]**.5 e**]j@ (eqn 3.18) 
V(f) = [real? + imaginary? ]**.5 e**j8 (eqn 3.19) 
iV(f) | = [real2 + imaginary2 ]**.5 "(eqn 3.20) 
since je**j@] = 1 

Iv(f)1 = ein Geeeniates 4 (eqn 3.21) 


2. Properties of the Fourier ITfeansforn 


Several properties of the Fourier transform are 
useful in the study of digital signals. They are 
represented here without proof and without a great deal of 
detail. 

a. Linearity 
if vi(t) <-> V1(f) and 
if v2(t) <-> V2(£) then 
Fiavi(t) + bv2(t) J] = avV4(f) + bv2(f£) 
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Figure 3.4 Amplitude Spectrum of Square Wave 
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Figure 3.5 Phase Spectrum of Square Wave 


b. Time-delay 
Fiv (t-toy ] = V(£) e**-juto 
Note that the amplitude spectrum of the delayed 
versicn is the same as the amplitude spectrum of the 
undelayed version. It is comforting to note that the 
Fourier transform of a pulse is the same tomorrow as it is 
today. 
c. Scale change 
Fiv(aty} = Wat vesa) 
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Frequency translation 
v (t)cos(2 pi fe t) <-> 1/2[V(f+fc) + V(f-fc) ] 
v(t) is any voltage 
cos(2 pi fe t) is a carrier wave of frequency fc 

Since understanding of this very important property 
of the Fourier transform is essential to the understanding 
of digital signal modulation it is expanded slightly here. 

As we have seen, the Fourier representation of a 
voltage pulse of amplitude A and duration T is the sinc 
function of amplitude AT as illustrated in Figure 3.6. 
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Pigure 3.6 Sguare Wave in Tiswe/Sinc Function in Frequency 


The translation is the product of v(t) and in this case 
cos(2 pi fc t). In the time domain this product only exists 
in the interval between -T/2 and T/2 as in Figure 3.7. The 
significance of this translation and its relationship to the 
bandwidth of tne voltage will be discussed further in the 
section deaiing with bandwidth. 
e. Differentiation 
d v(t) /dt <-> jw V(f£) 
A differentiator could be used as a clock for timing 


but would never be used in the presence of noise since the 
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oo Figure 3.7 Hultiplication and Translation Diagrams 


a result is exaggerated in the presence of high frequencies 
ik because w= 2 pi f. . 

ae f. Integration 

ae v(t) at <-> 1/w ¥(£) 

‘ere An integrator could be used to reduce the effects of 
: noise since w=‘*2 pi f is in the denominator tending to 
= deenphasize the presence of high frequency noise. 
Ss B. THE SAMPLING THEOREM 

- Essential to the understanding of digital communications 
a is the sampling theorem which was first introduced by 
of Nyquist in 1928 (Bef. 9], and later by Shannon in 15948 
ae (Bef. 10jJ. The sampling theorem states tnat any voltag2 can 
Pate, be uniguely represented by appropriately spaced sanpie 
ite values of the original voltage. “ore correctly stated, the 

ne Sampling theorem places limits on tne accuracy with wnick a 
aoe signal can be represented. 

to The implication is that an analog signal can be 
ae represented digitally or by a set of numbers, i.e., Saavdies. 


face A description of the sampling theorem follows. 
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Given any analog voltage, v(t) as in Figure 3.8, the 
sampling theorem says that the entire analog Signal is not 
required to accurately represent the voltage but only 
samples of it, call them vs(t). Figure 3.9 shows samples of 
a representative analog voltage. Samples can be taken of 
v(t) at every T seconds for a period of seconds. This can 
be accomplished by the use of a voltage clock, call it 
ve(t). The sampling can be viewed grapnically as a block 
diagram representing an analog voltage multiplier as shown 
in Figure 3.10. To be of further use in the understanding 
of digital communications we are interested in a frequency 
description of the sample voltage, vs(t). Note that the 
system which describes the obtaining of vs(t) involves a 
voltage multiplication. It was demonstrated in the 
proceeding section that voltage multiplication amounted to 
frequency translation, a property of the Fourier transforn. 


v(t) 
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Figure 3.8 Representative Analog Voltage “3 

First it is necessary to find the frequency description — 

of the clock, vco(t), Let vc(t) be a periodic square wave of a 

height 1 and duration das in Figure 3.11. Since vc(t) is ei 

periodic, it can be shown that’ the Fourier series Re 
representing vc(t) is given by equation 3.22. 
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Figure 3.9 Sasples of a Representative Voltage 
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Figure 3.10 Analog Voltage Multiplier 


Again it is noted that vs(t), the sample voltage, is the 
product of v(t), the original analog voltage, times vc(t), 


the clock voltage. In other words vs(t) = v(t) times a 
series cf cosine terms. 
co 
vc(t) = ao +} an cos(2 pi o fe t) (eqn 3.22) 
ned 


ao and an are left unevaluated 
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Figure 3.11 A Representative Voltage Clock 


00 
vs(t) = v(t){ao + ) an cos(2 pin fe t) ] (eqn 3.23) 
nai 


Now, assuming any form of the Fourier transform of v(t), 
as in Figure 3.12 and since vs{t) isa product and by 
utilizing the frequency translation property of the Fourier 
transform, the following is oktained as a representation of 
the frequency spectrum of vs(t). See Figure 3.13. 

Remember that the curve highlighted in the box in Figure 
3.13 is the Fourier represention of v(t), the original 
Signal. This original signal can now be recovered by 
filtering with an appropriate low pass filter of bandwidth 
equal to or greater than Be. This low pass filter wcecuid only 
permit the reception of that portion of the signal which 
represents the original voltage. 

The only question lett to resolve is how often to take a 
sample. Again referring to the diagram in Figure 3.13, it 
is noted that in order to prevent any overlap of successive 
translations (aliasing) 
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Figure 3.12 Fourier Transform of v(t) 
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Figure 3.13 Pourier Transform of vs (t) 


fo - B >> B 

fo >> 2B 
or the frequency of the clock, fc (the sampling rate = 1/T) 
must be strictly greater than 2 times the largest 
Significant Fourier frequency component present. Since 8 
will vary with different v(t), the sampling rate necessary 


to uniquely recover that v(t) also varies. 
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In summary on the sampling theorem, it can be said that 
any analog signal can be represented as certain sample 
values spaced the appropriate distance apart. The sampling 
theoren places limits on the accuracy of this 
representation. These sample values can then be transmitted 
from one position to another using analog to digital 
conversion and any one of a variety of modulation 
technigues. It has been shown that the Fourier frequency 
representation of the sample is not eguivalent to the 
Fourier frequency representation of the source voltage; 
however, the source voltage can be recovered at tae 
receiving erd by filtering. A block representation of the 
system is shown in Figure 3.14. 
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Figure 3.14 Analog Voltage Multiplier and Low Pass Filter 


C. THE AUTOCORRELATICN FUNCTION 


The above concept certainly seems simple enough. Ifa 
Signal is present and it is desired to transmit it froa 
point A to point B, all that has to be done is to take 


approfriately spaced samples of the signal, encode thea, 
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somehow transmit thea to point B, filter the received signal 


and the uniqueness of the original voltage has been 
reproduced. Unfortunately it isn't quite that easy. The 
reason it is not, is due to the presence of noise. The 
sources of noise will not be discussed here, however, the 
effects of noise in general terms and how a faint signal can 
be recovered in the presence of noise will be discussed. 

The time varying descriptor of voltage and the frequency 
varying descriptor of voltage have been introduced. Both 
are precise mathematical de scriptions of something (voltage) 
that is deterministic. That is, it can be described in 
mathegatical or graphical terms during any period of time. 
Such is not the case for noise corrupted voltages which are 
entirely randon. Therefore other descriptors of voltages 
must be employed in the presence of noise. ‘They are partial 
descriptors of voltage since a random signal cannot be 
described with precision. These partial descriptors are the 

1. Autocorrelation function and the 
2. Probability density function (p.d.f.) 

It is important here to note the difference between the 
source voltage at the receiver and the sample voltage. The 
source voltage at the receiver, vsr(t), is the digitally 
converted sample voltage, vs(t), modulated onto a carrier by 
one of the digital modulation techniques yet to be 
discussed. Understanding the autocorrelation function is 
the basis for understanding how a decision is made that 
vsr (t) is present at the receiver in the presence of noise. 
It shouid be remembered that the exact form of the carrier 
is kncewn at both the transmitter and receiver. 

First of all, an additive noise nodel is assumed where 


the signal at the receiver, vri(t), is equal to the signal 
which is being transmitted, vsr(t), (i.@e, the digital 
samples of v(t) modulated onto the carrier), plus randon 


noise, vn(t) 


+e te tp tt 


vr (t) 


vsr (t) + vn (t) (egn 3.24) 


ver (t) 
vsr(t) = signal voltage at the receiver 


voltage at the receiver 


vn(t) = noise voltage at the receiver 


If it is realized that most receivers operate on the 
principal of detection of DC voltage, the problem becomes 
one of rectification of the received signal and 
determination if the DC voltage, characteristic of the 
transmitted signal is present. 

A common type of rectification of an analog voltage 
involves squaring the input waveforn. If the given signal 
at the receiver is vr(t) as in equation 3.24 above, squaring 
the waveform introduces the square of not only the desired 
Signal cut also the sguare of the noise term. If instead of 
squaring vr(t) and introducing or at least not eliminating 
the neise, vr(t) is multiplied by vsr(t), the results in 
equation 3.25 are obtained. 


vo (t)evsr(t) = vsr*#2(t) + vsr(t) vn (t) (egn 3.25) 
By averaging, all that remains is vsr*¥*2(t) since the 
average of vsr(t)vn(t) is zere because vwn(t) is random and 


the average of any random voltage is zero. The DC component 
of vr(t)evsr(t) is vsr ° Any remaining AC component 
can be removed by a low pass filter. Graphically in block 
diagran form this receiver looks like Figure 3.15 
illustrated below. If vsr**2(t) > 0 then vsr(t) is present 
in the signal vr(t). If vr (t) is pure noise or some other 
Signal is present exclusively, then the output of the 
receiver will be 0. 
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Figure 3.15 Rectification through an AVM and LPF 


In practice, however, the exact waveform of vsr(t) may 
not be known or it may be a time delayed or distorted 
versicn of the original signal when it arrives at the 
receiver. If the distortion or delay is significant enough 
there will be little agreement or "“correlaticn" in the 
receiver. ; 

The solution is to multiply the received signal by a 
series of time delayed approximations of the transmitted 
signal. The signal is present whenever the output of the 


series of voltage multipliers or correlators exceeds a 


certain threshold approaching vsr**2(t). The concept is 
illustrated in Figure 3.16. 

For the concept illustrated in Figure 3.16 to work, the 
average value cf vr(t)evsr(t-p) must have a DC’ component. 


VC (t) evsr(t-p) = vsr(t)vsr(t-p) + (egn 3.26) 
vn (t) vsr (t-p) 


Vr (t) evsr (t—p) VSr(t) VS (t-p) + (eqn 3.27) 


vn(t) vsr (t-p) 
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Figure 3.16 Correlation Device 


But the average of vn(t}vsr(t-p) = 0 since vn(t) is random and 


the results of equation 3.28 are obtained. 
vr (t) evsr(t-p) = vsr(t) vsr (t—p) (eqn 3.28) 


What is of interest is the average value of 
vsr (t)vsr(t-p). The autocorrelation function (ACF) will be 
defined as that average value of vsr(t)vsr(t-p). The 
mathematical representation of the autocorrelation function 
for a continuous voltage is represented in eguation 3.29 
while the autocorrelation function for a puise voltage is 
presented as equation 3.30. 


a 
ACF = Rvv(p) = 27 f{ v(t) v(t-p) at (eqn 3.29) 
00 -T 
ACF = Cvv(p) = f v(t) v(t-p) at (eqn 3.30) 
-©o 


The autocorrelation function is a measure of the degree 
to which two identical signals which are corrupted in some 
manner are alike. The ACF of two signals which are 


identical, not distorted in any way and occurring at exactly 
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the same time, i.e, p = 0 has a maxinum value. On the 
other hand there will be very little correlation between two 
identical time variant signals when the time difference 
between then is great. In this case the autocorrelation 
function is near zero. 

A similar concept is used within the context of signal 
comparison. This concept is the crosscorrelation which is a 
measure of the degree to which 2 different signals are 
alike. When the crosscorrelation between voltage vsr and 
random noise voltage vn is 0 there is no relationship or 
correlation. 

Again let us assume that the signal at the receiver is a 
noise disrupted version: of the signal originally 


transmitted. 
vr (t) = vsr(t) + wn (t) (eqn 3.31) 
In the receiver , the time delayed version of the signal is 


applied to the incoming signal and the average is formed as 
before and shown in Figure 3.17. 
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Figure 3.17 Voltage Correlator 
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{ vsxr(t) + vwn(t) Je{ vsr(t-p) + vn(t-p) ] = (eqn 3.32) 
vsr (t)vsr(t-p) + vn (t)vsr (t-p) + 
vsr (t)vn(t-p) + vn (t)vn(t-p) = 


RvSI vsr(p) + Ryn vsr(p) * Rvsr va(p) + Rvn vn (p) 


u 
oO 


Rvn vsr(p) 


RvScr va(p) 0; because average vn = 0 


Therefore the result is simply Rvsr = vsr(p) + Rvn vn(p) 
or the autocorrelation function of the desired signal plus 
tne autocorrelation function of the noise. Assuming the 
shape of the autocorrelation function of both components 
(both the original signal and the noise which vary with p, 
time)is known, what is done in practice is to vary the value 
of p until the ratio cf Rvsr ver(p) / Rvn vn(p) is a maxinuo 
or the autocorrelation function of the signal is a maximum 
while the autocorrelation of the noise is minimum and tae 


signal is recovered. 


D. THE MATCHED FILTER 


The matched filter is a linear filter which has the 
characteristic response desired for optimum reception of the 
desired signal. In other words, we desire the response of 
the system to the linear filter to be in some way 
proportional to the autocorrelation function of the desired 
signal and in no way related to the autocorrelation function 
of the noise. This system could be illustrated as in Figure 
3.138. 
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= vsrlt)+ vn (+) a 


o& Rysr VSr (+) 
Figure 3.18 Matched Filter Block Representation 


a The question is what should the makeup or h(t), (i-e., 


= the response of the linear filter) be? The output of a 
a - linear system is a convolution so the output of the linear 
ss filter, h(t) to an input vsr (t) would be vsr(t) * h(t), 
: where * denotes convolution. 

In the previous Section it was stated that the 
a autocorrelation function, Rvsr vsr(p), a function of p, is ‘ 
Sy given by . 
me) yore) 
i RVSC vsr(p) = [vse (t)vsr(t-p) dt (eqn 3.33) 
act — 0o . 
oe By simply changing variables, equation 3.33 can be rewritten 
as a function of time, t. 
a co 
a Rvscr vsr(t) = [vst (ey vsx (p-t) dp (eqn 3.34) 
“.- — 0d 
> From convoluticn theory it is known that the response of a 
~ linear filter as above, vsr(t) * h(t) can be expressed 
= mathematically as in equation 3.35. 
y oo 
- vsr(t) * h(t) = [vsr(p) h(t-p) dp (eqn 3-35) 
a aay 
aia 46 
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As stated earlier, the desired output of the linear 
filter is the autocorrelation function Rvsr vsr(t). Note 
the similarities between the last two equations. If in 
equation 3.35 a simple substitution of h(t) = vsr(-t) is 
performed, the desired results are obtained. 

fore) 


vsr(t) * h(t) = [vsx(p) vsr (p-t) dp (eqn 3.36) 
—COo 


= RVSC vsr (t) 


Thecefore, for every signal that is desired to be 
recovered, the matched filter will do the job very nicely if 
the response of that filter, h(t) is equal to the inverse of 
the signal that is desired to be detected. 


47 


wp ele sires © 5 cage e 
oe on ge ss » wi o 
a Pe . 


a: 2 Pees 2 4%! 


IV. ANALOG TO DIGITAL CONVERSION 


Key to the understanding of digital communications is 
oa how the analog information is converted to digital 


ae, 4° 


information. As illustrated earlier, each analcg voltage or 
_ Signal can be represented by appropriately spaced sample 
; values. These sample values are still analog in that they 
can take on any value in the range of the original analog 
signal. What is desired is to change this infinite set of 
decimal numbers to a finite set of decimal numbers. This is 
: called analog to digital conversion. Common types of A to D 
ms conversion include Pulse Code Modulation (PCM), Differential 
a Pulse Code Modulation (DPCM), Delta Modulation (DM), Pulse 
‘ Amplitude Modulation (PAM), Pulse Duration Modulatior 
(PDM),and Pulse Position Modulation (PPM). PCM, DPCM, and 
DM will ke examined in this chapter because of their 
widespread usage and easy application to digital technology. 


A. PULSE CODE MODULATION (PCM) 


PCM is the most widely used A to D conversion technique. 

It involves assigning the sample value obtained in sampling 
‘ to one of several quantization levels within the range of 
the voltage sampled and then representing those quantization 
: levels by various binary code words. 
For example, the highest significant frequency component 
a in human speech is 3300 Hz. In order to capture accurately 
2 the signal produced in speech, one would require a sampling 
rate greater than or equal to 2B or 2(3300) = 6600 
samples/sec. The telephone company uses 8000 samples/sec. 
r B = 3300 Hz 
f 2 2B = 6600 samples/sec 
use £ = 8000 samples/sec to ensure no aliasing 


- 48 


PM ak RT a ee ee eee Oe RT Re es aF 


‘ats 

ot If a signal is sampled at the rate of 8000 ‘samples/sec, then 
bes in the transmission of those samples, assuming the samples 
a are transmitted immediately and not delayed, the 
ee transpission rate must be 

we Ts = 1/f 

a, Ts = 1/8000 samples/sec 

- Ts = 125 microsec/sanple 

be The number of quantization levels employed to represent 
i the various analog sample values is arbitrary. The greater 
~ the number of levels, the more accurate the reconstruction 


of the original signai but the more rapid the data 
transmission rate must be(in all cases there will be some 
& error present after recovery). Consider again the exaaple 
: of voice as illustrated in Figure 4.1. 
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Th 2125 nsec 
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Figure 4.1 Samples of Voltage 
2 
ce Assame it is desired to represent each sample value (decimal 
a numoer) with an 8 bit binary code word. Then the number of 
s quantization levels is given by eguation 4.1. 
2 The spacing between quantization levels is the range of 
Y voltages to be represented divided by the humber of 
A quautization levels. The analog to digital conversion takes 
: . 
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Figure 4.2 Analog Nature of Sarples 
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#guantization levels = 2**n (eqn 4.1) 


a for n = 8 bits/binary code word 
a #quantization levels = 2**8 
= 256 
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piace by determining into which quantization level the 
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respective sample values fall and assigning that sample 
value the binary number associated with that guantization 


oo 
See | 


RGN level. 

Now each sample value is assigned to an eight kit binary 
= code word. For voice traasmission it has teen shown that 
the Minimum transmission rate was Ts = 125 microsec/sanmple. 
The bit rate is then simply equal to the number of bits per 
sample times the inverse of the transmission rcate as 
illustrated in equations 4.2 and 4.3. 
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If voice was to be time multiplexed onto a_ single 


roe 
soe 


Channel, the bit rate for that channel would be the number 
of signals per channel times the bit rate for the most tiae 
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intensive signal as illustrated in equation 4.4, 
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bit rate = (n Lits/sample) (1 sample/T sec) (eqn 4.2) 
bit rate = (8 bits/sample) (1 sample/125 microsec) 

= 64 kbps 

bit duration = 1/bit rate (eqn 4. 3) 


= 1/64 kbps = 15.6 microsec/bit 


(*signals/channel) (bit rate/signal) = (egn 4.4) 
bit rate/channel 


In other words, if 24 voice signals are to be 


multiplexed onto a single channel, the bit rate must be 


bit rate/channel 2 (eqn 4.5) 
(24 signals/channel) (64 kbps/signal) 
= 1536 Mbps 


AS can be seen, the data transmission rate increases 
Significantly, necessitating better and better hardware and 


a wider and wider bandwidth. Advantages should be readily 


apparent for A to D conversion technigues or modulation 


technigues that reduce the bit rate required. 


Recovery of the analog signal is by a process called 


Digital to Analog Conversion. All that will be said about D 


to A conversion is that it is the inverse process of A to D 


conversion and that a slight error is always introduced 


during the process due to the discrete nature of the 


quantization process during A to D and D to A conversion. 


DIFFERENTIAL PULSE CODE MODULATION (DPCM) 


In CPCM, what is converted to binary numbers is not the 


quantization level (decimal number) which each sample value 
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is represented by but the successive differences between 
quantization levels. The idea is that the range of maximum 
difference values will be smaller than the range of actual 
sample values. It is therefore possible to represent that 
range of delta values with fewer bits/binary word and 
therefore fewer quantization levels and ultimately a lower 
bit rate. [Ref. 11] 

In actual practice, in a DPCM conversion technique, it 
is a statistical estimate of each successive sSampie value 
which is subtracted from the actual value that is ccnverted 
to binary code words. The result is the same in that the 
range of amplitudes is reduced and therefore fewer bits/word 
are reguired to represent the sample thereby lowering the 
data transmission rate required to transmit the signal. 
(Ref. 11) 


C. DELTA BODULATION (DA) 


Delta modulation is a form of DPCM where successive 
quantization levels of the output differ by only 1 bit. 
That is to say that any successive quantization level can be 
represented by varying only one bit of successive output 
binary code words. In other words a type of gray code is 
employed. This A to D technigue is implemented through the 
use of a DM coder or linear delta modulator. This DM coder 
approximates a given input signal with a series of linear 


segments of uniform slope. A comparison is made between the 
value of this approximation andthe input signal at each 
sample increment. The Sign of this difference value is what 
is encoded andis used to increment the DM coder in the 
direction of the input signal. By using the differential 
sign value of the input signal and the incremented 
approximation from the DM coder the linear approximation 
from the linear delta modulator is said to "track" the input 
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signal. { Ref. 8] Slope overlcad of this type of modulation 
technigue occurs when the slope of the incoming signal 
exceeds the ability of the DM system to follow *he source at 
the sampling rate being utilized. [Ref. 11]. 
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V- DIGITAL SIGNAL HODGLATION LECHEIOUES 


Digital signal modulation techniques are the methods of 
encoding information for transmission utilizing digital 
technology. Factors affecting digital modulation include, 
but are not limited to, the physics of the method, hardware 
requirements, bandwidth considerations, power requirements, 
data transmission rates and error probabilities. It is 
these aspects which will be explored in the following 


secticns. 


A. DIGITAL HODULATION FORMATS 


Modulation is the technique by which the characteristics 
or one waveform (called the carrier) are varied or nodiried 
by the characteristics of another (called the source). The 
carrier waveform of interest is the sinusoid. It should be 
obvious taat the attributes of a particular sinusoid that 
differentiate it from every other sinusoid are its 
amplitude, phase and frequency. It follows that the 
characteristics of the carrier waveform to be varied by the 
source are its amplitude, phase, frequency or any 


combination of the three. This gives rise to the broad 


general formats of Amplitude Shift Keying (ASK), Phase Shift 
Keying (PSK) and Frequency Shift Keying (FSK). All other 
digital modulation technigues are variations of or 
combinations of these basic formats. Other factors involved 
in the description of the digital modulation technique keiag 
employed include the number of bits being encoded at one 
time, the employment of error correction technigues and the 


baseband (source) waveforn. 
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When each bit of the baseband waveform is individually 
encoded by any of the techniques previously mentioned (ASK, 
PSK, FSK), the technique being utilized is referred to as 
binary encoding. When more than 1 bit of the source code is 
~ modulated onto the carrier at one time it is called block 
encoding. Block encoding allows for one of m = 2**k 
waveforms where k = ausber of bits. 

This paper will net go deeply into signal detection or 
- demodulation techniques, however, a basic understanding of 
what is involved is necessary and again further defines the 
digital nodulation format being employed. In general terms, 
Signal detection is referred to as either coherent or 
noncoherent detection of the transmitted signal. Coherent 
detection, perhaps the easiest to understand, is when all 
possible waveforms of the modulated carrier waveform are 
available at the receiver and the waveform at the receiver 
is in phase with the transmitted carrier. Noncoherent 
detection is involved when the receiver does not have 
knowledge of the phase of the transmitted information and 
we one of a number of phase estimation techniques must be 
employed for Signal recovery. 


B. HARDWARE 


Althougn significant advances have been made in recent 
i years to improve the quality of hardware associated with 
Satellite communications, most components are not what would 
be considered “off the shelf items", The hardware 
components most commonly referred to with regard to digital 
communications include the sampler, encoder, nodulator, 
oa multiplexer and transmitter. There are variations of the 
| above hardware requirements necessary for certain types of 
~ digital formatting, however, those exceptions will be 


addressed separately when the individual modulation 
technigues are examined. 


Figure 5.1 Basic Hardware Component Block Diagraa 


he 1. Sampler 

é The sampler is the device or component which samples 
é or extracts those single characteristic values of a 
ae haturally occurring analog signal which are ultimately 

encoded into a digital word by the encoder. The frequency 
at which this sampler must operate was derived in the 
discussion of the Sampling Theoren. It was determined that 


a the sampler must operate at a frequency greater than 2 times 
Pes the highest significant Fourier freguency component of the 
source voltage. Modern solid state devices capable of 


taking thousands of samples/sec are available at modest 


costs. 
2. Encoder 


Often referred to as the A to D converter (analog to 
digital converter), the encoder transforms the samples of 
the analog signal derived from the sampler into a diyital 
format through one of the analog to digital conversion 
technigunes described in the chapter on A to D conversion. 
These digital bits can be stored for later use, coded, 
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delayed or used immediately either individually or in groups 
to modify one of the characteristic qualities of the carrier 
waveforn. When sample values of the analog source signal 


are ccnverted into digital bits of information, they are 
simply that, digital bits of information. The carrier can 
only be modified to represent the source information by 
interaction with another voltage. Therefore the value of 
the digital bit (binary), either 0 or 1, is used to generate 
a basektand waveform or voltage which does the actual 
modulation of the carrier waveform. 


a. Common Baseband Waveforms 


It should be obvious that since it is desired to 
represent binary code words with a representative taseband 
waveform what is required is two levels of voltage. There 
are two basic logic schemes for representation of the 
baseband waveforn. They are bipolar or unipolar logic. 
Bipolar logic involves representing the 0's or 1's of the 
binary codeword as either +V or -V as illustrated in Figure 


5.26 
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Figure 5.2 Bipolar Logic 
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The 0 or the 1 of the binary code word can be represented 
either as +V¥, -¥ or -V¥, +V respectively. This is a matter 
of convention but must be clearly understood in the various 
component designs in order to ensure compatibility between 
parts of the systen. 

Unipolar logic utilizes a voltage to represent 
either of the possible binary digits and the absence of 
voltage to represent the other. Two common conventions of 
unipolar logic are represented in Figures 5.3 and 5.4. 


+ 
Le | 


Figure 5.3 Unipolar Positive Logic 


Again, which form of unipolar logic might be employed in the 
A to D converter is a matter of convention. 

Variations of these two basic encoding formats 
have the advantages of ease of generation and improved 
decoding and clock recoverability. Two variations commonly 
used in satellite communications are the Non Return to Zero 
(VEZ) and the Manchester waveforms. The NRZ waveforms is 
simply the voltage representation which corresfonds to the 
stream of bits represented in the logic scheme chosen. 
There is no transition as long as the same bit is present. 
Choosing bipolar logic as an example, a NRZ waveform can be 
described schematically as in Figure 5.5. 
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Figure 5.4 Unipolar Negative Logic 
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Figure 5.5 Non Return to Zero Waveforr 


The Manchester waveform of the same source 
voltage is represented with a transition at the midpoint of 
the bit duration from either +¥ to -V or -vV to +#V. This 
form cf coding has the advantage of clock synchronization in 
all cases of digital encoding. Schenatically, the 
Manchester code can be represented as in Figure 5.6. 
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Figure 5.6 Manchester Waveforn 
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odulator 


The modulator in a digital communications system is 
more commonly referred to as a mnoden 
(modulator/demodulator). This device does the actual 
transformation of the digital information into a waveéforn 
that can be transmitted from one point to another. The 
modulation as mentioned earlier can ‘modify either the 
amplitude, phase or frequency of the carrier wave. In 
Simple terms, the digital information derived from the 
analog source is mapped into the carrier wave. It is this 
mapping that determines the power and bandwidth 
Characteristics of the moden. 


4. Multiplexer 


Multiplexing is a common method for increasing the 
utility of a given communications channel. Variations in 
multiplexing technigues give rise to the multiple access 


techniques employed in satellite communications. 
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a. Prequency-Division Multiplexing (FDM) 


Frequency-Division Multiplexring is a technique 
whereky the capacity of a communications channel is 
increased by adding one signal to another. These signals 
occupy discrete nonoverlapping frequency bands. A specific 
signal is recovered by use of a band pass filter for the 
frequency band desired. [Ref. 7] 


be. Time-Division Multiplexing (TDM) 


Tine-Division Multiplexing is a technique used 
to increase the capacity of a given communication channel by 
adding signals in time. That is, specific signals or 
portions of signals are allocated specific time slots or 
portions of the carrier wave. These time allocations cannot 
overlap and the signal is recovered through proper 
synchronization of the recovery hardware with the 
multiplexer. [Ref. 7] 


c. Code-Division Multiplexing (CDM) 


Code-Division Multiplexing is a technique for 
increasirg the capacity of a given communications channel 
through the assignment of a characteristic code to the 
digital signal before it is multiplexed in the time or 
frequency domain. Since the code used in multiplexing is 
known at the demultiplexer, recovery of the digital source 
code could be accomplished through the use of a nmatched 
filter or correlator. [{Ref. 7] 


C. BANDWIDTH 


For purposes of discussion and for uniformity, when 
examining the various digital sigaal modulation techniques 
of interest, the bandwidth will be defined as the highest 
Significant Fourier frequency component of the signal in 
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question. For digital communications, it can be shown that 
by the above definition, the bandwidth of the baseband 
waveform is 1/Tb, where Tb is the bit duration. For 
exasple: 
1. Voice sampled at the rate of 8000 samples/sec 
2. Each sample represented by an 8 bit ccde word 
3. Rb = bit rate = (8000 samples/sec) (8 bits/sample) 
Rb 64 kbps 
4. Tb bit duration = 1/Rb 
Tb 15.6 microsec 
The Fcurier transform of the baseband waveform cf duration 
15.6 microsec is represented in Figure 5.7. 
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Figure 5.7 Pourier Transform of Baseband Waveform 
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It should te readily apparent that as the bit rate of the 
digital baseband waveform increases, so does the bandwidth. 
It is desireable to keep the bandwidth of the signal toa 
minimum due to limitations on availability of the 
electromagnetic spectrum. In addition, interference with 
Signals which occupy adjacent frequency bands is reduced and 
propogation limitations of the. medium and hardware 
limitations such as self-inductance and capacitance are 


minimized. 


D. SBECIFIC TECHHIQUES 


Now with some of the basics of digital signals already 
introduced it is possible to investigate some of the more 
significant modulation formats employed in digital satellite 
communications. The organization of this section will be 
around the three general types of modulation, i.e., phase, 
amplitude and frequency. Within each broad category, 
several techniques will be described under the subareas of 
binary encoding, where the information is modulated onto the 
carcier bit by bit, and block encoding, where groups of bits 
do the modulation. Variations of these areas will be 
pointed out. The characteristics of the modulation 
technique will be introduced including, where practical, an 


analytical description, phasor representation, error 
probability, spectral efficiency, advantages and 
disadvantages. 


Phase shift keying is a generic form of signal 
modulation where the baseband waveform is used to modify or 
change the phase of the carrier. Phase, as stated earlier 
is one of the characteristics which distinguish one sinusoid 
from another. In general, phase shift keying offers the 
advantages of power and bandwidth efficiency. 


a. Binary Phase Shift Keying (BPSK) 


As the name suggests, Binary Phase Shift Keying 
(BPSK) results in a mcdulated carrier waveform consisting of 
two phases or the same sinusoid. The baseband waveform is 
used to change the phase of the carrier bit by bit from one 
phase to the other. In the case of bipolar logic, the bit 
stream consists of +1. The general waveform of the carrier 
can be represented as a cosine wave of amplitude A, angular 


63 


oe. we 


frequency w and initial phase angle delta as in equation 
5.1. ‘ 


vc(t) = A cos(wt + delta) (eqn 5.1) 


Since two phases are represented in vc(t), it is 
custosary to let then differ ty pi radians. This can be 
represented by the scdificaticn of equation 5.1 to account 
for a pi phase shift depending on whether the bit, v(t) to 
be sodulated is either a +1 or -1 as Shown in equation 5.2. 


vco(t) = A cos(wt + delta + pi/2(1-v(t))) (egn 5.2) 


Equation 5.2 can be expanded to resuit in a simplified fora 
showing that the modulated carrier waveform is simply a 
product of the baseband waveform and the carrier. See 
equation 5.3 where v(t) represents the baseband waveforn. 


ve(t) = A v(t) cos(wt + delta) (eqn 5.3) 


A phasor diagram can be constructed to represent 
the modulated carrier waveform in BPSK and is shown in 
Figure 5.8. As can be seen by the phasor diagram in Figure 
5.8, the phase of the modulated waveform defends on the 
value of the tktaseband waveform and differs by pi depending 
on whether the value of the modulated bit is +1. 

Remenbering the section on the properties of the 
Pourier transform, it is possible to analyze the spectral 
and power efficiency of BPSK. Recall BPSK can be created by 
multiplication. Multiplication of two voltages is covered 
by the frequency translation property of the Fourier 
transfora. Frequency translation results in double the 
bandwidth or spectral requirement of the translated wavefora 
and half the power. In the case of a NRZ baseband wavefora, 
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Figure 5.8 Phasor Diagram of BPSK 


the bandwidth of the non-translated pulse, determined from 
the largest significant Fourier frequency component is 1/Tb, 
where Tb is the bit duraticn. The bandwidth of the 
translated voltage (BPSK carrier) is 2/Tb. Tke power in a 
sinuscid is A2/2 and none of this power capacity is lost in 
BPSK since only the phase cf the original sinusoid is 
changed. For a Manchester baseband waveforn, the 
non-translated signal has bandwidth 2/Tb and the translated 
wave has bandwidth 4/Tb. 

The advantage of BPSK is that it is relatively 
efficient in bandwidth utilization and power capacity. It 
does have the disadvantage of fairly large sidelobe 
components which contribute to interference with adjacent 
frequencies. This usually necessitates some type of 
filtering before the signal can be transmitted. 


b. Differential Binary Pnaase Shift Keying (DBPSX) 


The recovery of a Binary Phase Shift Keyed 
signal requires the use of a phase coherent reference signal 
as described in the previous sections on autocorrelation and 
matched filters. This may not always be practical or 
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& possitle soa technigue called Differential Binary Phase 
a Shift Keying has be developed. DBPSK uses the same carrier 
= type as BPSK and offers many of the same advantages and 
Ps disadvantages. As stated, the primary difference is the 
ne lack of necessity for a phase coherent reference. 
Bie The baseband waveform, v(t), for DBPSK is 
ae generated by comparing the phase of the next bit to be 
ie transmitted to that of the previously transmitted bit. For 
‘Se example, if a +1 is encoded onto the carrier as v(t), it 
Oe represents a particular phase cf the carrier. If the next 
i bit to be sent is also a +1, there is no phase change, i.e., 
a the bit remains the same over two succesSive bit durations. 
= If, however, the second bit to be sent is a -1, this 
> represents a pi phase change and is decoded as a bit change 
“ over successive bit durations. This can ke expressed 
a mathematically as in equation 5.4 where bk represents the 
= bit which will be encoded on the carrier and bk-1 represents 
i the previous bit already encoded and ak represents the 
“ present bit. 
aan bk = bk-1 e ak (eqn 5.4) 
. As can be seen when bk = -1, there is a phase 
pas change and the bit changes from one bit to the next 
LO depending on the logic type being employed. If bk = #1, 
this represents the phase remaining constant over two 
successive bit duraticns, i.e., no bit change. See Figure 
= 5.9 [Ref. 8]. Recovery is effected Ly correlation of the 
me received bit to the previously received bit. 
ce Differential Binary Phase Shift Keying (DBPSK) 
- differs from Differentially Encoded Binary Phase Shift 
a Keying (DEBPSK) in the recovery step ouly. DPPSK takes 
oe advantage of the modulation of phase shift in the recovery 
Se process while DEBPSK still utilizes a phase coherent carrier 
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Figure 5.9 Differential Binary Phase Shift Keying 


sigral in demodulation. Both cases use the same scheme for 


differential encoding. 


c. Orthogonal Binary Phase Shift Keying (Orthogonal 
BPSK) 


In Binary Phase Shift Keying there is always a 
probability, with the interjection of noise on the 
transmitted signal that an individual bit will be decoded 


incorrectly from that which was actually transmitted. In 


some environments, this probability is so high that special 
encoding technigues must be employed in order to reduce the 
probability that a bit or a tinary code word will te daecoded 
in error. One such technigue is Orthogonal BPSK. 
Orthogonal BPSK is a type of group encoding scheme where 
qroups or blocks of bits from the original signal are 
assigned redundant bits according to a predetermined 
assignment routine. This new expanded sequence of bits is 
thei transmitted as before in BPSK. 
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For example, assume that it has been determined 
that the respective values of the original signal will be 
represented by a k bit binary code word. It has been 
demonstrated that there exists 2**k binary code words, k 
bits long in this type of arrangement. With Orthogonal BPSK ; 
there also exists 2**k orthogonal sequences of bits (called 
channel symbols, bauds or Chips) [Ref. 6], which represent 


the k binary code words. However, the number of bits or 
chips in the sequence is increased to n= 2**k, thereby 
providing the desired redundancy. Each of the 2**k 


sequences of n chips iS constructed to be orthogonal to 
every other sequence in order to assure maximum separation. 
See Figure 5.10 [ Ref. 6]. 
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Figure 5.10 Orthogonal Binary Phase Shift Keying 


Decoding is aiso done by blocks through a bank 
or 2**k correlators. Since the decoding is dore in blocks, 


the prokability that the entire transmitted code word will 


of 
si 
oda 


be incorrectly decoded is reduced even though individual 


aaa 
ene 


a fad te 


bits are incorrectly decoded. This is the major advantage 
of Orthogonal BPSK. Figure 5.11 [Ref. 6], shows the 
reduction in error probability with Orthogonal BPSK 
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Figure 5.11 Probability of Error for Orthogonal BPSK 


The major disadvantage to Orthogonal 
that it results in a reduced data bit rate at 


BPSK is 
a given 


transmission rate or the necessity for a higher transmission 


rate if the same bit rate is to be maintained. For example, 


suppose the bit rate of a PSK Signal is R bits/sec. 
of the k bits inthe PSK signal are representea 
chips as is the case in Orthogonal BPSK modulation, 
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If each 
by 2**k 


then the 


bit rate must be (2**k/k)R if the same data rate is to be 
maintained. Additionally, Since the necessary bit rate 


increases, the bit duration decreases. It was shown in the 
secticn on bandwidth that as bit duration decreases as the 
result of a higher bit rate, an increase in bandwidth is the 
result. Therefore, the tradeoff in decreased transmission 
error comes at the expense of bandwidth and Lit rate. 
(Ref. 6] 


d. Quadrature Phase Shift Keying (QPSK) 


Quadrature Phase Shift Keying derives its name 
fron a four-phase modulation of the carrier waveform. These 
modulations are achieved by simultaneously modulating the 
inputs from two bit streams onto a single carrier. The two 
bit sequences could be from two separate sources or 
successive bits from a single source. To take advantage of 
the orthogonality of the sine and cosine functions, QPSK can 
be viewed as in equation 5.5 as the sum of a BPSK modulated 
sine and cosine. 


ve(t) = A vi(t) cos(wt + delta) + (eqn 5.5) 
A vaz(t) sin(wt + delta) 


The signal can also be represented as in 
equation 5.6 where the phase angle, theta, of the modulated 


carrier is shown in eguation 5.7. 
vco(t) = 2**.5 A cos(wt + theta + delta) (eqn 5.6) 
theta = arctan (v1i(t)/vz(t)) (egn 5.7) 
The phases possible as a result of substitution of the 


possitle bits in bifolar logic in equation 5.7 are +45 
degrees and +135 degrees. See Figure 5.12. 
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Figure 5.12 Sodulaticn and Phases of QPSK 


Since two bits are being encoded at one time for 
transmission, the bit rate in QPSK is twice that of BPSK. 
The BESK bit rate is equal to the QPSK symbol rate, (symbol 
rate being the time the quadrature bits exist on the 
carrier). It is the symbol rate in QPSK which determines 
the bandwidth and since symbcl rate is equal to bit rate in 
BPSK, the bandwidth of the two modulation techniques is 
ideatical. The great advantage of QPSK lies in the 
bandwidth utilization. Two Lits of information are 
transmitted thereby effectively doubling the bit rate at no 
additional cost in bandwidth. A phasor representation of 
QPSK is illustrated in Figure 5.13. Each quadrature 
component of the modulated signal contains half the power of 
the tctal carrier or (A2/2)/2 = A2/4. 


e. Offset Quadrature Phase Shift Keying (OQPSK) 


*: a aa ve a i 7 


In Quadrature Phase Shift Keying, Loth bits of is 
the two baseband waveforms, which are to be modulated onto 1 
the carrier, change at the beginning of the respective = 
symbol time. This allows the phase of the carrier to change I 


sa 


a 


up to 180 degrees fcr each symbol as illustrated in the 
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Figure 5.13 Phasor Diagram of QPSK 


phasor diagram in Figure 5.13. The phase change of the 
modulated signal can be limited to 90 degrees through , a. | 
modulation technique called Offset Quadrature Phase Shift 
Keying. This is accomplished by delaying the application of 
the second baseband channel for Ts/2 seconds, where Ts is 
symbol duration as illustrated in Figure 5.14. Modulation 


is then accomplished as before in QPSK. 
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Figure 5.14 Offset Quadrature Phase Shift Keying 
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bes 

es The result is a QPSK modulated signal in which 
Ke the maximum phase shift between successive bits is 90 
; degrees. The bandwidth and power spectra of an OQPSK 
cat : modulated signal are the same as that of a QPSK modulated 
aa signal. The advantage of OQFSK is in the limit which is 
3 placed cn phase shift during encoding which simplifies the 
sl encoding hardwa.-e. Additionally OQPSK has spectral and 
ae interference advantages during decoding [{Ref. 6]. 
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an f. Multiple Phase Shift Keying (MPSK) 

ms Another form of digital modulation is known as 
oe Multiple or M-ary Phase Shift Keying (MPSK). Once again, as 
es in Orthogonal Binary Phase Shift Keying, k data bits are 
pa transmitted. The k bits represent the word length and there 
Ss are m = 2**k different binary woris, k bits long, for this 
oe type of modulation. In MPSK, the binary code word is used 
» to vary the phase of the carrier. There are therefore m = 
a 2**k different phases in MPSK. These sequences of binary 
me digits need not be orthogonal as in Orthogonal Binary Phase 
os 
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Shift Keying. The length of the transmitted symbol is only 
k bits long as compared to 2%**k bits for Orthogonal BPSK. 
The phases therefore are not all orthogonal leading to 
greater difficulty in decoding and a higher probability that 
an individual code word will be incorrectly decoded, 
especially as m becomes large. 

As in other types of Phase Shift Keying, MPSK 
can be represented as in equation 5.8 where theta = phase = 
(2pi/s) i; i=0,1,-.-m-1. 


vc(t) = A cos (wt + theta + delta) (eqn 5.8) 


An advantage of MPSK is that as long as the 
symbol transmission rate does not increase, the carrier 
bandwidth does not increase even as the size of the binary 
code word increases. The disadvantage, as stated earlier, 
has tc dec with error rates in decoding. Decoding requires a 
phase coherent reference of considerable Stability 
especially as m increases and respective phases become 
closer together. Practical limits for this type. of 
modulation is m = 8 or a3 bit binary code word due to the 
complexity of the enccding aud decoding hardware. 


2. Applitude Shift Keying (ASK) 


Amplitude Shift Keying, as the name implies is 
modulation of the carrier through variation of its amplitude 
due tc variations in the baseband waveform. The information 
in the baseband waveform is thereby imparted to the carrier 
through this modulaticn. Since the amplitude of the carrier 
varies between successive symbols or binary code words, the 
power also varies. Recovery is by comparison of the 
transmitted power during each symbol to the possible ao = 
2**k power levels. 
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ae* Multiple Amplitude Shift Keying (MASK) 


MASK is a type of group encoding where k bits 


(amp 


* are combined into a single waveform for transmission and 
: subseguent recovery. Since each symbol or binary code word 
- contains k bits, there are m = 2**k different symbols anda 
consequently mn different amplitudes of the carrier possible 
during the symbol duration time. An assignment. scheme which 
maps the binary code words into the m different amplitudes 
would be to simply space them Ai = A/mi volts apart, where A 
is the maxinum amplitude and mi represents the decinal 
equivalent of the binary code word from i = 1 to 2**k. 

The analytic form of a MASK waveform is 
represented in equaticn 5.9 where all the variables are the 
: same a presented earlier except for Ai which is equal to one 
A of i different amplitudes depending on the bit sequence to 
be transmitted. 


vce(t) = Ai cos(wt + delta) . (eqn 5.9) 


MASK is not very popular in satellite 
communications since it depends on a carrier of very stable 


amplitude, This is not generally practical under actual 
. coniitions. However, MASK does have the spectral advantage 
sa of a constant bandwidth of 2/Ts, TS equal to symbol 


duration, even as the binary code word increases in length. 
be. Quaarature Amplitude Shift Keying (QASK) 


Quadrature Amplitude Shift Keying (QASK) is a 
a hybrid digital signal modulaticn technique. It represents 
. separate amplitude modulation of the quadrature components 
of a ccmmon carrier. In that sense it is both ASK and PSK. 
Implementation of QASK involves simultaneous application of 
_ M-ary Amplitude Shift Keying to the quadratures. See 


Ca ee A 


equation 5.10, where Al and A2 are derived according to an 
assignment scheme as in MASK. 


vco(t) = Ail cos(wt + delta) + (eqn 5.10) 
A2 sin(wt + delta) 


By this method, the effective bit rate is 
doubled over that of ordinary MASK. In other words, QASK 
results in the same data rate as modulating 2k bits onto the 


carrier at the same time or during the same symbol duration 
with MASK. The advantage to QASK lies in the fact that this 
increase in bit rate comes at no further increase in 


bandwidth. The disadvantages associated with MASK are still 


present with an additional increase in the complexity of the 


decoding equipment. Decoding of QASK must be done in two 


steps. First the signal is recovered in phase through a 


*. phase-coherent correlator and then each amplitude modulated 
signal is recovered as before in MASK by comparison of the 
power levels of the received signal. 


ms 3. Frequency Shift Keying (FSK) 
a 


Frequency Shift Keying (FSK), is the generic tera 


used to describe a number of digital signal modulation 


techniques which cause variations in the carrier frequency 


: by interaction with the baseband waveforn. Decoding is 


generally through measurement of the frequency of the 


received signal. 


a a. Mimimum Shift Keying or Fast Frequency Shift 


Keyiag 


Power requirement for the transmission or 


retransmission of a satellite signal is extremely important. 


< It is desireable to limit power required to accurately 


Se oe - 
: < ak 


transBit a piece of information to the minimum possible 
consistent with allowable error rates. For this reason, 
detection of FSK signals is limited to coherent methods as 


power required increases significantly for non-coherent 
methods. 
ae MSK or FFSK are identical and represent a 
i modulation technique known aS continuous phase FSK. They 
get their names from the fact that "E£ast" indicates more 
bits per second can be transmitted in a given bandwidth than 
ordinary BPSK and “*minimum" refers to the minimum nodulation 
index for which orthogonal signalling occurs [Ref. 7}. 
Analytically, FFSK can be expressed as in 
equation 5.11 and equation 5.12. 


ve (t) A cos((we + deltaw) t) (eqn 5.11) 


vc (t) 
A sin(tdeltaw t) sin(we t) 


A cos(+deltaw t) cos(wc t) - (eqn 5.12) 


As is noted, FFSK can be envisioned as_ the 
summation cf separately modulated in-phase aod quadrature 
components of the carrier by the baseband waveforn. The 
technique is Similar to that used in OQPSK. The frequency 
variation in the carrier waveform is between woe + deltaw and 
wo ~ deltavw. Maximum separation in phase of these two 
freguency components occurs at pi as noted in equation 5.13, 
where Tb represents the bit or symbol duration. 


(w2-w1)Tb = pl (egn 5.13) 


This can be converted to the modulation index 
menticned above simply by converting to frequency as in 
equation 5.174. 
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h = (f2-£1)Tb = .5 (eqn 5.14) 


For FFSK, the frequency deviation or separation 
between frequencies of the transmitted carrier is exactly 
1/2Tb and the deviation from the carrier frequency is 1/4Tb. 
FFSK can therefore be rewritten from equation 5.12 as in 
equation 5.15 by substitution for deltaw = 2pi(1/4Tb). 


vc(t) = A cos(tpi t/2Tb) cos(2pi fc t) - (eqn 5.15) 
A sin(tpi t/2Tb) sin(2pi fc t) 


The error rate in this modulation technique is 
identical to that of BPSK. Frequency Shift Keying can be 
implemented to incorporate any given modulation index by 
simple calculation of the frequency separation with equation 
5-14 and substitution into eguation 5.15. For M-ary 
Prequency Shift Keying, the separation between frequencies 
must ke at least 1/Ts, where Ts is the symbol duration, in 
order to avoid carrier energy from one frequency being 
incorrectly interpreted as carrier energy from another 
frequency. This results in a modulation index of 1. Figure 
5.16 represents the phasor diagram of an FFSK modulated 


Signal. 
4. Quadrature Partial Response Signalling (QPRS) 


Quadrature Fartial Response Signalling (QPKS) 
represents a specific type of the general class of digital 
signal modulation techniques called Partial Response 
Signalling. Partial Response Signalling is a method in 
which a controlled amount of intersymbol interference (ISI) 


is allowed during the encoding process. Previously any 
overlap in successive signals resulted in aliasing and an 
increased probability of error in decoding. The idea was 
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Figure 5. 16 Phasor Diagram of FFSK or MSK 


first introduced by Lender in 1963 [Ref. 12]. His concept 
was that in knowing and controlling the amount of 
interference allowed during the encoding process, 
compensations can then be made for the ISI at the receiver. 
Allowing for a limited amount of ISI makes it possibile to 


transmit at rates eyual to the Nyquist rate, something that 


is not possible with many other systems [Ref. 13]. 

QPRS is implemented, aS with previous quadrature 
systems, thrrough simultan: cus modulation of the guadrature 
components of a common carrier. However, this time the 
modulation is accomplished with a PRS ‘systen. The 
modulation of the respective quadrature components 
represents the impulse response of one of a number of linear 
filters to the bits cf the baseband waveform. The linear 
filter employed depends on the class of the Partiai Response 
Signalling system being used. There are 5 classes of linear 
filters commonly used in PRS and they are illustrated 
analytically in Table 1 and graphically in Figure 5.17. 

Output of a QPRS system can be represented as in 
equation 5.160, where an and bn represent the nth bit to be 
modulated and h(t - oTs) represents the contribution of the 
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be Figure 5.17 QPRS Linear Filter Impulse Responses 
impulse response of the linear filter at time t for bit n. 
a {Ref. 14] 
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ve(t) = )an h(t - n Ts) cos(wt + delta) (eqn 5.16) 
Ar. 
kn h(tec n Ts) sin(wt + delta) 


n2x0O 


QPRS has the advantage of rapid transmission rates 
with no increase in bandwidth and excellent error handiing 
performance. The cost of such performance improvement comes 
in the form of a higher required signal to noise (SNR) ratio 
when compared to other binary systems. 
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Due to the nature of digital signal modulation, it is 


possikle to simulate systems through the use of the digital 
computer. This chapter is a description of the coaputer 
Simulation included in the Appendix of this thesis. 


A. GENERAL DESCRIPTION OF THE PROGRAM 


The program was constructed using top-down design and 
modular programming. FORTRAN was selected as the 
programming language due to its capanilities in the area of 
numerical operations andthe availabiiity of mathematical 
routines to be used as modules in the progran. Testing was 
accomplished on each module as the program was developed. 
The program was designed to be used interactively in order 
to allow for instructional use. Although the program was 
made user friendly in that error checking is accoygplishei on 
user inputs, care must be taken to insure instructions are 
followed correctly. 

The. program consists of twenty modules, sixteen of which 
were written by the author, three which were taken from the 
Double Frecision International Mathematics and Statistics 
Library (IMSLDP) and one which was taken from a NON-IMSL 
divrary which resides on the IBM 3033 iocated at the Naval 
Postgraduate School. The IBM 3033 was used exclusively for 
development and testing. 

The development was accomplished using the WATFIV 
compiler; however, the program was written to operate with 


the VS FORTRAN compiler as well. Testing and operation was 
done using the VS FORTRAN compiler. Figure 6.1 is a block 
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diagram of the relationshif of the modules for a 
representative modulation subroutine. 


Be. MAIN CONTROL AODULE 


MAIN operates as the control module for the entire 
computer simulation and accomplishes limited output. It 
introduces the programs and allows the user to input various 
parameters of the digital signal modulation technique to be 
Simulated and various other control functions. These inputs 
include: 

1. Modulation technique 

2. Class of QPRS system (when appropriate) 
3. Digital logic scheme 

4. Baud or symbol rate 

5. Bits/binary code word (when appropriate) 
6. Carrier frequency 

7. \Number of samples to be generated 

8. .Carrier max amplitude 

9. Initial phase angle 

16. Use or random number generator or no 

11. Seed for RNG (when appropriate) 

12. Number of repetitions of the simulation 

Once the user has input the desired characteristics of 
the modulation to be simulated, MAIN calls the appropriate 
module to begin the actual calculation. The values 
necessary to calculate the time series of the signal are 
passed to the subroutines as parameters. MAIN calls the 
reguired subroutine the number of times the user specifies 
as repetitions. Each call to a subroutine produces the 
Discrete Fourier Transform and amplitude spectrum of the 
signal. It is these respective amplitude spectrums which 
produce the statistics in the final output after MAIN calls 
the statistics subroutine the final time. 
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Figure 6.1 Representative Block Diagraaz 
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C. SUBROUTINE BPSK 


This description of SUBROUTINE BPSK, a module that 
Simulates Binary Phase Shift Keying, will serve as a guide 
on the construction and operation of subsequent modules 
involved in the actual signal modulation. A complete 
description including its interaction with other program 
subroutines will be included. Since other signal modulation 
subroutines interact in the same manner within the frogran, 
subsequent descriptions will concentrate on differences to 
this kasic module. 

SUBROUTINE BPSK begins with variable declarations as 
with all FORTRAN programs. All variables are declared to be 
either type integer or double precision. All calculations 
in the program are carried out in double precision accuracy. 

Variable initialization is next. Significant variable 
initialized at this point are TIME, the time of the first 
sample which is set to 0. STEP, the delta time between 
samples is set to the normalized -Nyquist sampling rate. 
OMEGA, the angular frequency and DELTA, the initial phase 
angle are computed in radians. NBAUD, a variable which 
keeps track of the number of symbols which have been 
nodulated, is set to 1 or the first symbol is being 
modulated... NBAUD remains at its present value until total 
elapsed time exceed 1 symbol duration or BAUDD, the 
normalized baud rate. 

Subsequently, an initial call is made to SUBROUTINE 
STREAM to receive the value of the first bit to be 
modulated. The modulation is then carried out according to 
equation 6.1. 

The values of each sample are assigned to an array for 
storage along with the corresponding time increment. 
Modulation continues with the first drawn bit until 1 bit or 
symbol duration is exceeded when STREAM is called to draw 
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oe 
oa vc(t) = A v(t) cos(wt + delta) (eqn 6.1) | 
oe 
aes 
‘x A = amplitude 
= v(t) = bit to be modulated 
Rot: w = angular frequency 
Es t = tine 
, delta = initial phase angle 
Sa 
aa another bit. This entire process is repeated until the 

number of samples specified by the user in MAIN is reached. 
is At this point, only on the first repetition of the 
aoe Simulation, SUBROUTINE PLOT is called which gives the user 
ee the opportunity to plot the time series of the simulated 
oe Signal. Upon return from PLOT the subroutine calls the IMSL 

: routine FFTRC to generate the Discrete Fourier Transform of 
: : the time series produced. This is accomplished on each 
Sis repetition of the program. 

; Again on the first repetition only, information on the 
ae humber of the principle harmonic of the FFT is displayed. 
Sent Plot is called again to give the user the opportunity to 
oe plot the amplitude spectrun,. The values of the FFT and the 
a amplitude spectrum are computed on each successive 
pots repetition of the program and the amplitude spectrum is 
oy added to the statistics being accumulated by a call to 
By SUBROUTINE STATS. 

Oa 

ie D. SUBROUTINE DBPSK 

ee SUBROUTINE DBPSK simulated Differential Binary Phase 
oe Shift Keying and is essentially the same as SUBROUTINE BPSK. 
aes The difference lies in the fact that what is encoded during 
ae the uodulation process is the product of the random bit 
oe drawn and the value of the Lit previously modulated. A 
Deu 
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Ss reference bit equal to +1 is used to determine the first bit 
r to be nodulated. In this manner, modulation of +1 means no 
change occurs between successive bits and the phase of the 
carrier remains the same. Modulation of a -1 indicates a 
4 change of bits has occurred and is indicated by a change of 
k phase of the carrier between successive bit durations. The 


rest of the module remains unchanged. 


E. SUBROUTINE OBPSK 


Orthogonal Binary Phase Shift Keying is accomplished in 
this module. Again SUBROUTINE OBPSK is constructed 
essentially the same as SUBROUTINE BPSK. In Orthogonal BPSK 
the only bit streams that are modulated are n bits long 

% representing each binary code word, where n = 2¥**k (k = 
number of bits per word is limited to 6 bits in MAIN). Sach 


sequence of n bits is orthogonal to every other allowed 
sequence. 

The way this is accomplished in this module is first to 
draw a series of k random bits from SUBROUTINE STREAM. This 
series of k 1*"s and O's is’ then changed to its decimal 
equivalent from 0 to 2**k-1 and saved for later use. 

The program continues with a call to SUBROUTINE ORTHO 
which generates an n X nm orthogonal matrix of +1's and -1'"s 
through the use of the Hadamard matrix and the Kronecker 
product [Ref. 15]. The Hadamard matrix is a 2 X 2 
orthogonal matrix of 1*s and -1's. The Kronecker product is 
the matrix which is formed when a matrix is expanded to 
‘ twice its original size by replacing each element of the 
: Hadamard matrix by the product of the Hadamard matrix and 
: the original matrix. The value of the decimal equivalent to 
the k random bits is used to identiry the row of the matrix 
and the respective columns represent the value of the binary 
digit which is modulated onto the carrier by BPSK. This 
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Sequence continues until either one bit duration is 
exceeded, when another orthogonal bit is received from ORTHO 
or when all the orthogonal bits in a row are modulated. Then 
another sequence of k random bits is converted to its 
decimal equivalent and a new row and column of the 
orthogonal matrix is identified. Execution continues in 
this manner until all required samples have been produced. 
It should be noted that the bit rate for Orthogonal BPSK is 
n times the baud rate. This significantly increases the 
Signal bandwidth reguirement and decreases the time 
increment between successive samples in the Simulation. The 
rest of the module is the same as those previously 
presented. 


F. SUBRCUTINE QPSK 


Quadrature Phase Shift Keying, Simulated in SUBROUTINE 
QPSK, is accomplished through simultaneous BPSK modulation 
of the guadrature components of the carrier. Two random 
bits are drawn from successive calls of SUBROUTINE STREAM 
and each bit is in turn modulated onto the carrier 
components and the sum formed. This process is repeated at 
time intervals equal to the normalized Nyquist sampling rate 
until TIME exceeds one symbol duration, BAUDD. At this time 
two new random bits are drawn and the modulation continues 
until all required samples are produced. The rest of the 
modules is the same as those before. Eguation 6.2 
represents the analytical expression used to sinulate 


G. SUBROUTINE OQPSK 


SUBROUTINE OQPSK is essentially the same as SUBROUTINE 
QPSK with these minor deviations. The time cr the first 
sample is artificially initialized as TIME = .5(BAUDD) or 
1/2 the symbol duration. In this manner, the two randon 
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vc(t) = A vi(t) cos(wt + delta) + ; (eqn 6.2) 
A v2(t) sin (wt + delta) 


vi (t) 
v2 (t) 


in-phase bit 


quadrature Lit 


bits first being modulated are both known to have existed on 
the carrier for 1/2 the symbol duration. The first of thest 
two bits only remains on the carrier until TIME = 1(BAUDD) 
when a third random bit is drawn and modulation begins with 
this bit. The second random bit is allowed to remain until 
TIME = 1.5(BAUDD) when a fourth random bit is drawn to 
replace it and so on. In this manner the two bits being 
modulated are offset in the time they change by 1/2 the 
symbol duration. This necessitates keeping track of the 
number of bits modulated on both the in-phase and quadrature 
component of the carrier. The remaining portion of the 
module is the same as SUBROUTINE QPSK. | 


H. SUBROUTINE MPSK 


This module simulates M-ary Phase Shift Keying. This is 
accomplished by changing the phase of the carrier to any one 
of n =2**k phases where k is the number of bits in the 
binary code word. The max length of the binary code word is 
limited to 10 inthe MAIN prograa. This modulation is 
accomplished by use of equation 6.3 to simulate the MPSK 
systen. 

The program first draws k random bits and converts then 
to the decimal equivalent. At this point modulation begins 
at TIME = 0 according to equation 6.3. Modulation continues 
until one symbol duration is exceeded when a new set of k 
random bits is drawn, conversion to decimal takes place and 
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ve(t) = A cos{wt + (2pi m)/n + deltaj (eqn 6. 3) 


decimal equivalent of binary code word 


2**k; k = bits in binary code word 


nodulation continues or until all desired samples have been 
produced. The programs continues as in previous subroutines. 


I. SUBROUTINE MASK 


The principles cf operation of SUBROUTINE MASK is 
Similar to SUBROUTINE MPSK which also involves a group 
encoding system. SUBROUTINE MASK allows for M-ary Amplitude 
Shift Keying. The analytic expression used in the 
Simulaticn is provided in equation 6.4 Ai is one of 2**k 
equally spaced amplitudes. 


vc (t) = Ai cos(wt + delta) (eqn 6.4) 
Once again a set of k random bits is drawn from 
SUBROUTINE STREAM. The decimal equivalent of the k bits is 
computed and used to adjust the amplitude Ai accordirg to 


the assignment routine expressed in equation 6.5. 


Ai = A/sm (eqn 6.5) 


max amplitude 
m= 1tons;s n = 2**k 


This process is repeated at each symbol duration until 
all desired samples have been produced. The rest of the 


module remains unchanged. 
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J. SUBROUTINE QASK 


SUBROUTINE QASK combines the elements of the previously 
described amplitude modulation technique with the now 
familiar quadrature sodulation technique. Two separate 
streams of random numbers are generated and converted to 
their decimal forn. These nunbers are used to produce two 
separate amplitudes with the same assignment scheme used in 
MASK. These amp itudes modulate the quadrature components 
of a common carricc and the output is formed by the sua of 
the guadratures. Modulation continues in the above manner 
until a symbol duraticn elapses when two new amplitudes are 
calculated. Processing terminates when all required samples 
are ccaputed. 


K. SUBROUTINE ASK 


Minimum Shift Keying is frequency shift keying in which 
the modulation index is 1/2. The modulation index is 
represented in equaticn 6.6 


h = .5 


(deltaf) Tb (eqn 6.6) 


deltaf = frequency separation 
Tb = bit duration 


The frequency deviation trem the carrier is .5(deltaf) = 
1/4Tb. Converting this freguency deviation to radians 
yields pi/2Tb. This angular frequency deviation from the 
central carrier is either +pi/1Tb depending on the value of 
the bit to -ve modulated. 

The simulation in SUBROUTINE MSK is accomplished by 
drawing a random bit and generating the signal according to 
equation 6.7. 
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vc(t) = A cos{ (w +pi/2Tb)t + delta] {egn 6.7) 


Once again modulation continues until one bit duration 


. is elapsed when another bit is drawn. Modulation stops when 
~ all reguired samples are _ produced. The rest of the 
re subroutine remains unchanged. 


L. SUBROUTINE MFSK 


This subroutine modulates a signal simulating M-ary 
Frequency Shift Keying. The separation between adjacent 
freyuencies is established as 1/Ts, where Ts is the symbol 
duration. This amounts to a moduiation index equal to 1. 
Initially the ertire range of frequency deviation from the 
carrier is calculated as (n-1)TS. .The mean frequency 
deviation is then calculated. In other words the range of 
frequencies of the modulated signal is the frequency of the 
Re central carrier +1/2 the magnitude of the entire range of 
. frequency deviation. The Minimum frequency is used as the 
base for computing the frequency to be used to modulate the 


r Signal during each respective symbol duration. 

Pt: At this point in the subroutine a series of k randon 
2 bits (k<10) 1S drawn and converted to decimal. This decimal 
J number is multiplied by the frequency separation and added 
ae to the Minimum frequency and converted to radiars. The 
a analytical expressicn of the signal produced with this 
w subroutine is illustrated in eguation 6.8. 

: vc(t) = A cos{(w + 2pi mf)t + delta) (eqn 6.8) 
‘s mx = modulation freguency 
= 


This sodulation continues until a symbo* duration 


“ elapses when a new frequency deviation or mf is calculated 
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or mecdulation stops do to completion of all required 
sanples. 


Hf. SUBROUTINE QPRS 


This mcdule simulates 5 different classes of Quadrature 
Partial Response Signalling systems. The technique employed 
involves the summation, at ¢€ach time increment, of all 
responses to the linear filter representing the class. This 
is accomplished for all of the bits being modulated during 
the duration of the simulation. 

Initally, two arrays of random bits are generated 
representing those bits which would be modulated onto the 
quadrature components of the carrier during the length of 
time the signal is to be simulated. The the impulse 


_ response to each of these n bits is calculated for the time 


of the respective sample. These impulse responses represent 
the modulation for the bit in question onto the respective 
porticns of the carrier and the sum is formed according to 
equation 6.9. 


fos) 
BS (t) = Da h(t - nis) cos(wt + delta) + (eqn 6.9) 
n=-Go 
Tes h(t - ofs) sin (wt + delta) 
nz- 


h(t - nTs) = impulse response of nth bit at time t 


The time of the first sample is initialized at 6(BAUDD) 
in order to ensure that when the first sample is taken, 
contributions from bits modulated at time 0 are also summed 
do to the overlapping nature of the QPR signals. The value 
of the last sample is also formed by the sum of the 
quadrature components existing at 6(BAUDD) after it is 
produced. The rest of the program operation remains the 
same aS previous modules. 
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BS. SUBROUTINE STREAS 


This subroutine interacts with the IMSL random number 
generator SUBROUTINE GGUBFS to produce a randoa nunber 
between 0 and 1 and make assignment on the basis of the 
value of this random number to random bits according to the 
digital logic scheme specified in the MAIN program. It also 
enables the user to manually insert binary digits if that 
option was specified previously in MAIN. 


O. SUBROUTINE ORTHO 


CRTHO produces an n Xn orthogonal matrix from the 
Hadamard matrix by forming successive Kronecker products 
{Ref. 15]. It also selects the appropriate row and column 
of the orthogonal bit to be sodulated and passes this bit 
back to SUBROUTINE OEBPSK. Foint of entry to SUBROUTINE 
ORTHO is controlled by a flag set in OBPSK. 


P. SUBROUTINE PLOT 


PLOT is an interactive module that allows the user to 
determine whether or not a graph of the output of program 
calculations is to be produced. The user is also able to 
selectively vary certain parameters associated with the 
plot. PLOT calls SUBROUTINF UTPLOT which actually produces 
the graph and performs the output. UTPLOT is a NON-IMSL 
library routine from the Naval Postgraduate School. 


Q- SUBROUTINE STATS 


The Statistical calculations associated with the 
ampiitude spectrums of the various functions generated on 
successive rcepetitions of the program are computed in 
SUBROCTINE STATS. Upon completion of each repetition, each 


module calls STATS where the sus, sum of the squares, sum of 


: the cubes and sum of the quartics of each element of the 


ae amplitude spectrum of the signal produced are computed and 
we stored. When the last repetition of the program is 
ae complete, MAIN instructs STATS to compute and output the 
a final statistics, i.e., mean, vactiance, skewness and 
oe kurtosis of each component of tke amplitude spectrun 


according to the estimations contained in [Ref. 16]. 
ye Additionally the mean variance and variance of the variance 
ie is calculated and output. Point of entrance to this 
a subroutine is controlled by flags set inthe individual 

modulation subroutines or in MAIN. 
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VII. STATISTICS OF IHE FAST FOURIER TRANSFORMS 


A. DESCRIPTION OF THE PROBLEM 


As stated in Chapter II, if it can be shown that the 
statistics of the FFT can be somehow linked to a particular 
digital signal modulation technigue, then the nodulation 
technique can be identified on the basis of those statistics 
alone. Although it is beyond the scope of this thesis to 
actually derive those relationships if they exist, it seen 
prudent to attempt to determine if there is statistical 
differences between the components of the FFT's as they are 
derived in the enclosed computer simulations. 

The statistic chosen to do the hypothesis testing is the 
F-test since the true mean and variance of the distribution 
need not be known [Ref. 17]. In addition, the information 
necessary to calculate the F-statistic is readily 
accumulated in SUBROUTINE STATS. Calculation of the 
F-statistic necessitated the writing of a small computer 
program to be used for that purpose once the output from the 
Main frogram was generated and reformatted. 


B. ABHALYSIS-OF-VARI ANCE 


If there were no differences attributable to the 
modulation technique, then it would be reasonable to assume 
that for a certain set of characteristics, the components of 
successive FFT's would be the same. In other words, if a 
signal was modulated by BPSK and a statistically significant 
number of FFT's were generated of the signal, then the 
statistics of those FFT's would be assumed to be related. 
On the other hand, if it can be shown that the statistics 
are not related (i.e€., not from the same distribution), then 
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‘el significance can be placed in the variation among modulation 
@ techniques. 

an 1. Zhe F=Distribution 

«<2 

~! The assumptions necessary to use the F-distributiona 
» are: 


1. Normal distribution 
2. Random samples 


‘ 3. Independent samples 

me These assumptions are not too difficult to intuitively 
- accept in the model that will be proposed. Once again it 
nS wonld be expected that the FFI's of successive identically 
: modulated signals to be related. For the simulations 
a conducted as part of this test, they were all generated 
: using bits from a random number generator for which it can 


. be shown that each sequence of bits passes statistical tests 
a for randomness and independence. The assumption of 
‘a normality could also be argued on the basis of the central 
limit theorem. 


h , The value of the F-statistic with which comparison 
5 will be made is 1.51 for 14 degrees of freedcm in the 
3 numerator, an infinite number of degrees of freedom in the 
denominator and a 90% confidence region. How these values 
ag were olktained will be detailed under the design of the 
. experiment. 

§ 2- Design of the Experiment 


The experimental data used in the computation of the 
Statistics was chosen to minimize the random contributions 
of variables other than the modulation technigue. The 
computer simulations included in this thesis were used to 
generate the statistics and a separate program also included 
in the appendix was used to calculate the F-statistic. The 
variables in the experiment include: | 


a 7) Ag eT sty te 


se a 
ot 2 ot 


i 


al 


97 


‘os io a. - - “ - aoa + oe oe. ee ee er ee ee ee nF ~ ‘gs een ue 7 . o - e 7 . —~* ax o * ? at $4 8 o: 
ee ee re ee ee oa” eft et at ee go es gee ee ge cea a Os . 

: oat A ee en ee Sar ae es Ye 
PE CPC TLR rR ens Pata Pret TE Sp otigs Goh at en aeee See ate te 
oI - ‘a 2 A 5 x . ase. & Z bd - ¥ a id Z 


f 


ee hd 
OD Edt 


1. Modulation technique 

2. Logic typé 

3. Baud rate 

4. Bits per binary code word 

5. Carrier frequency 

6. Carrier amplitude 

7. Initial phase angle 

8. Number of samples generated 

9. Time between samples 

10. Seed for randos number generator 

11. Number of repetitions or trials 

Fifteen di fferent modulation techniques were 
compared. In all cases bipolar logic was simulated, the 
baud rate was held constant at 1200 baud ard the maxinua 
carrier amplitude was established at 1 volt. Additionaily 
the phase angle of all simulations was 0 degrees and the 
seed for the random number generator was 1 thereby ensuring 
the same random sequence of bits. The number of bits per 
binary code word did vary among the modulation techniques. 
When the modulation technique was M-ary, the bits per code 
word was always 3 so it is possible to infer that bits per 
code word is a function of the modulation technique. The 
time between samples was the normalized Nyquist sampling 
rate, Tais did vary between modulation techniques but was 
necessary to derive an accurate FFT. Carrier frequency also 
varied tetween Simulations but was always twice the lowest 
carrier frequency recommended in the computer simulation. 
In most cases this was 2400 Hz, or twice the baud rate. 
Finally each siaulation was repeated 100 times and the 
statistics of the FFT's based on those 100 repetitions. The 
Number of samples produced was always 64 so there are 33 
components of the FFT. 
What this amounts to can be illustrated through an 

Analysis-of-Variance table as shown in Table 2 [Ref. 17]. 
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TABLE 2 
ANALYSIS-OF-VARIANCE TABLE 


n Observations in Each of r Groups 


Meee er cease tn ct pe net came ae cee: ED cee nes ee dened thew nema ead) 


In the case of the simulation described above, this amounts 
to 100 observations (repetitions or trials) from each of 15 
groups (modulation technigues) for each of the 33 components 
of the FFT. The degrees of freedom in the numerator is 


equal to (15-1) or 14 woile the degrees of freedcm in the 


denominator is equal to 15(100-1) or 1485. Table values for 
the F-distribution use infinity as the degree of freedom tor 
values greater than 120. 

An F-test was also performed to determine if there 
is a relationship among the mean variance, mean skewness or 
mean kurtosis of the 15 different modulation techniques. 
Again the assumption is made that these statistics would be 
normally distributed from modulation techniques which had 
the same or statistically similar components of their FFT's. 
In this case there were again 15 modulation techniques to 
compare (14 degrees of freedom) but only 33 observations 
(the variance, skewness or kurtosis of the respective 
elements of the FFT). This number of observations yields 


15(33-1) = 480 degrees of freedom in the denominator. 
99 
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3. Hypothesis and Hypothesis Testing 


The hypothesis posed for the model is that the means 
of the individual components of the FFT's from the 15 
modulation techniques are from the same distribution. 
Therefore, the respective means must be equal. Also tested 
is that the mean variance, mean skewness and mean kurtosis 
of all components of the FFT's of each nodulation technique 


, 
‘a 
rn 


are equal. These hypothesis are illustrated in equations 
7.1 through equation 7.4. Each hypothesis is tested and 


AS ase ae 
ha! 2 in 


compared to the F-distribution selecting a rejection region 
or level of significance of .9. This equates to an 
F-statistic of 1.51 for 14 degrees of freedom in the 
numerator and an infinite number of degrees of freedom in 
the dencnminator. 


4. Result 


The means of the 33 ccmponents of the FFT for each 


of the 15 modulation techniques were compared using the 
F-distribution. The results are shown in Table 3. Table 4 
shows the F-statistics associated with the comparison of the 
mean of the variance, mean of the skewness and mean of the 
kurtosis for the reSpective grours. In- addition a complete 
summation of the results of the F-test are included in the 
appendix. 


Cc. CCHCLUSION 


The results of the F-test indicate that the aull 
hypothesis uust be rejected at the .9 significance level and 
the alternate hypothesis accepted that the means are not 
egual for FFT components 9-14. In addition the F-statistic 
for the comparison of mean variance, mean skewness and nean 
kurtosis all fall in the rejection area. 
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Ho: XBAR11 = XBAR21 = © @ e = XBARij (eqn 7.1) 


H1: XBAR11 # XBAR21 # © @ e # XBARij 


PP e Pe 


i = ith modulation technique 
= jth component of the FFT 


Hos MVAR1 = MVAR2 = © @ e 


MVARIS (eqn 7.2) 


H13 MVAR1 # MVAR2 # e @ @ # MVAR 15 


HO: MSKEW1 = MSKEW2 = e e © = MSKEW15 (eqn 7.3) 
H1: MSKEW1 # MSKEW2 # e © e # MSKEW15 

Hos MKUR1 = MKUR2 = e @e @e = MKUK15 (eqn 7.4) 
H1: MKUK1 # MKUR2 # e © @ # MKUR15 


This indicates that there is a difference between the 
Statistics of of the PFT's of the respective modulation 
technigue,. Since all the parameters used in the generation 
of these statistics were held essentially constant among the 
trials, it can be assumed that these differences are due to 
the modulation technique emrfloyed. It has not been 
established yet what those differences may be. This is an 
area where future research and study is warranted. 

The results of this test would seem to point to those 
components of the FFT which offer the best opportunity to 
develcp those relationships or differences. FFT components 
9-14 have okvious differences; however, FFT components 20, 
24 and 25 also have large fF-statistics but do not fall in 
the rejection area. These components may also be 
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TABLE 3 
FP-STATISTICS OF FFT COMPONENTS 


T 


FFT AMONG WITHIN OTAL 
CCMPCNENT GROUP SUM GROUP SUM OF SQUARES F-STAT 


1 0.121 17.69 17.81 0.727 
2 0.112 18.48 18.59 0.642 
3 0.108 20.65 20.76 0.557 | 
| 4 0.135 25.40 25.54 0.562 
5 0.172 32.74 32.91 0.558 
. 6 0.254 41.72 41.98 0.647 | 
S H 7 0.422 50.10 50.52 0.893 | 
l 8 0.962 82.83 83.79 1.231 
: 9 7.591 390.2 397.8 2.063 | 
: 10 49.09 1681.0 1730.1 3.097 | 
1 53.01 1732.1 1785.1 3.246 
| 12 61. 56 2071.7 2133.3 3.152 
: 13 45.73 1832.8 1878.6 2.646 
_ [. 14 3.864 250.5 254.4 1.636 
§ | 15 0.539 98.13 98.67 0.582 | 
¥ | 16 0.366 110.1 110.5 0.353 
§ { 7 0.504 151.0 151.5 0.354 | 
. 18 0.943 140.3 141.2 0.713 
: | 19 1.299 158.5 159.8 0.370 | 
: 20 2.671 200.7 203.4 1.412 
{ 21 1.853 226.3 228.2 0.868 | 
. | 25 2.277 23221 234.3 1.041 
: { 23 1.712 197.6 199.3 0.919 | 
: 24 2.512 181.9 184.4 1.476 | 
3 { 25 2.144 171.3 173.5 1.327 
4 | 26 1.257 143.7 145.0 0.928 
x | 27 1.233 113.9 115.1 1.149 
: e Soe ee | 
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TABLE 3 
F-STATISTICS OF FFT COMPONENTS (con't) 

| 28 0.858 115.8 116.7 0.785 

29 0.547 63.72 64.27 0.911 
| 30 0.356 55.77 56.13 0.677 | 
\ 31 0.183 41.41 51.59 0.469 

| 32 0.051 13.79 13.84 0.389 
L 33 0.031 10.92 10.96 9.394 | 
a _ 
| TABLE 4 
| F-STATISTICS OF THE MEAN VARIANCES, | 
MEAN SKEWNESS AND MEAN K3RTOSIS | 
{ { 
| AMONG WITHIN TOTAL | 
| GROUP SUM GROUP SUM OF SQUARZS F-STAT | 
VARIANCE 3.144 E5 3.755 56 4.070 E6 2.370 | 
| SKEWNESS 5.307 £10 7.031 E11 7.562 E11 2.588 
| KURTOSIS 8.280 E11 °1.262 £13 1.344 E13 2.250 
a = each 


interesting to examine, In addition, it should be noted 
that the statistics of the FFT associated with each 
respective modulation technique also display some strixing 
differences. These statistics may prove in some way to bea 


fingerprint of the modulation techniques themselves. 
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APPENDIX E 

RESULTS OF F-TEST 
F-STATISTICS FOR THE 1 COMPONENT 
SUS OF SUMS = 0.2091487882294994D 02 
SUM OF SUMS**2 = 0.1810655469453404D 02 
SUMS OF SQUARE OF SUMS = 0.4128097510679048D 02 
AMCNG GROUP SUM OF SQUARES = 0. 1211883136154515D C60 
WITHIN GROUP SOM OF SQUARES = 0.1769374494346614D 02 
TOTAL SUM OF SQUARES = 0. 1781493325708159D 02 
DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000C000000D Ow 


F-STATISTIC = 0.72650632143&7600D 00 


F-STATISTICS FOR THE 2 COMPONENT 

SUM OF SUMS = 0.2305792717805020D 02 

SUM OF SUMS#*2 = 0.1894701234015291D 02 

SUS OF SQUARE OF SUMS = 0.4663825126655468D 02 

AMCNG GROUP SUM OF SQUARES = 0.1119371755000360D 00 

WITHIN GROUP SUM CF SQUARES = 0. 18480602982748736D 02 

TOTAL SUM OF SQUARES = 0.1859256700298739D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.642474646050 14737D 00 


F-STATISTICS FOR THE 3 COMPONENT 

SUM OF SUMS = 0.2693734535827635D 02 

SUM OF SUMS**2 = 0.2124674921974186D 02 

SUM GF SQUARE OF SUMS = 0.5921566381358867D 02 

AMCNG GROUP SUM OF SQUARES = 0. 1084095881685183D 00 

WITHIN GROUP SUM OF SQUARES = 0.2065459258160598D 02 

TOTAL SUM OF SQUARES = 0.2076300216977449D 02 

DEGREES OF FREEDOM IN NUMERATOE = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000000000uUD QO4 
F-STATISTIC = 0.5567362242775778D 00 
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P-STATISTICS FOR THE 4 COMPONENT 

SUA OF SUMS = 0.3175688631627065D 02 

SUM OF SUMS**2 = 0.2621104542300214D 02 

SUS OF SQUARE OF SUMS = 0.8070513093726266D 02 

AMONG GROUP SUM OF SQUARES = 0. 1347180903696015D 00 

WITHIN GROUP SUM OF SQUARES = 0.2540399411362951D 02 

TOTAL SUM OF SQUARES = 0.2553871220399911D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 


F-STATISTIC = 0.5624997485041869D 00 

F-STATISTICS FOR THE 5 COMPONENT 

SUM OF SUMS = 0.3669104349531116D 02 

SUM CF SUMS**2 = 0.3380633740317695D 02 

SUM OF SQUARE OF SUMS = 0.1069682473356158D 03 

AMCNG GROUP SUM OF SQUARES = 0.1721940248396148D 00 
aITHIN GROUP SUM OF SQUARES = 0.3273665492982079D 02 
TOTAL SUM OF SQUARES = 0.3290884895466041D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.5579331866788867D 00 


F-STATISTICS FOR THE 6 COMPONENT 

SUM OF SUMS = 0.4160097878948431D 02 

SUM OF SUMS**z = 0.431323€642020094D 02 

SUB CF SQUARE OF SUMS = 0.1408215358007360D 03 

AMCNG GROUF SUM OF SQUARES = 0.254%4544005119440D 00 

WITHIN GROUP SUM OF SQUARES = 0.4172415106219358D 02 

TCTAL SUM OF SQUARES = 0.4197860546270552D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.6468757561623442D 00 
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- SUM OF SQUARE OF SUMS = 0.3070018365343552D 03 


F-STATISTICS FOR THE 7 COMPONENT 

SUS CF SUMS = 0.4810387092348111D 02 

SUM OF SUMS**2 = 0.5206409381611744D 02 

SUM OF SQUARE OF SUMS = 0.1964653659903027D 03 

AMONG GROUP SUM OF SQUARES = 0.4219987280210724D 00 

WITHIN GROUP SUM OF SQUARES = 0.5009944015621441D 02 

TOTAL SUM OF SQUARES = 0.5052143888423548D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGKEES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.8934632362547178D 00 


F-STATISTICS FOR THE 8 COMPONENT 
SUM OF SUMS = 0.5623255941290002D 02 
SUM CF SUMS**2 = 0.8590285622098203D 02 


AMCNG GROUP SUM OF SQUARES = 0.9619512065933320D 00 

WITHIN GROUP SUM OF SQUARES = 0.8283283785563848D 02 

TOTAL SUM OF SQUARES = 0.8379478906223181D 02 

DEGREES OF FREEDOM IN NUMEKATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISZIC = 0.1231824736914027D 01 ! 


F-STATISTICS FOR THE 9 COMPONENT 

SUM OF SUMS = 0.9931605295569187D 02 

SUM OF SUMS**2 = 0.4043531056175881D 03 

SUB OF SQUARE OF SUMS = 0.1416642417283150D 04 

AMCNG GROUP SUM OF SQUARES = 0.7590638589699638D 01 

WITHIN GROUP SUM OF SQUARES = 0.3901866814447566D 03 

TOTAL SUM OF SQUARES = 0.3977773200344562D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000G000000D 04 
F-STATISTIC = 0.2063499133280512D 01 
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F-STATISTICS FOR THE 10 COMPONENT 

SUS OF SUMS = 0.2029942628781924D 03 

SUS OF SUMS**2 = 0.1757564757819924D 04 

SUM OF SQUARE OF SUMS = 0.7655692250451960D 04 

AMCNG GROUP SUM OF SQUARES = 0.4908580866354581D 02 

WITHIN GROUP SUM OF SQUARES = 0.1681007835315404D OW 

TOTAL SUM OF SQUARES = 0.1730093643978950D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM Il DENOMINATOR = 0.1485000000000000D 04 
F~STATISTIC = 0.3057309684192641D 01 


F-STATISTICS FOR THE 11 COMPONENT 

SUM OF SUMS = 0. 2004773302247246D 03 

SUM CF SUES**2 = 0.1811875077561617D 04 

SUM OF SQUARE OF SUMS = 0.7980439863108266D 04 

AMCNG GROUP SUM OF SQUARES = 0.5301029200839379D 02 

WITHIN GROUP SUM OF SQUARES = 0.1732070678930534D 04 

TOTAL SUM OF SQUARES = 0.1785080970938928D 04 

DEGREES OF FREEDOMS IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.3246332537532912D 01 


F-STATISTICS FOR THE 12 COMPONENT 

SUM OF SUMS = 0.2210834652773109D 03 

SUE CF SUMS**2 = 0.2165859765218944D O04 

SUM OF SQUARE OF SUMS = 0.9414279932820438D 04 

AMCNG GROUP SUM OF SQUARES = 0.6155753358218842D 02 

WITHIN GROUP SUM OF SQUARES = 0.2071716965890740D 04 

TOTAL SUM OF SQUARES = 0.213327449 9472928D 04 

DEGKEES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D Ou 
F-STATISTIC = 0.3151731454585565D 01 
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F~STATISTICS FOR THE 13 COMPONENT 
SUM OF SUMS = 0.1979583434471791D 03 

SO# OF SUMS**2 = 0.1904701075347160D 04 

SUM OF SQUARE OF SUMS = 0.7185446672840862D 04 

AMCNG GROUP SUM OF SQUARES = 0.4572946290150773D 02 

WITHIN GROUP SUM OF SQUARES = 0. 1832846608618751D 04 

TOTAL SUM OF SQUARES = 0.1878576071520259D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F~STATISTIC = 0.26464786714598604D 01 


F~STATISTICS FOR THE 14 COMPONENT 

SUM CF SUMS = 0292277362726 14619) 02 

SUM OF SUMS**2 = 0.2600827156886254D 03 

SUM OF SQUARE OF SUMS = 0.9540726819081072D 03 

AMCNG GROUP SUM OF SQUARES = 0.3863985704619236D 01 

WITHIN GROOP SUM OF SQUARES = 0.250541992869 544 4D 03 

TOTAL SUM OF SQUARES = 0.2544059785741636D 03 

DEGREES OF FREEDOA IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000C000000D 04 
F-STATISTIC = 0.1635887377498235D 01 


F<STAIISTICS FOR THE 15 COMPONENT 

SUM OF SUMS = 0.6880926336210145D 02 

SUM CF SUMS**2 = 0.1018300362534556D 03 

SUM OF SQUARE OF SUMS = 0.3695062863415098D 03 

AMCNG GROUP SUM OF SQUARES = 0.5385863804584070D 00 

WITHIN GnhOUP SUM OF SQUARES = 0.9813497339004056D 02 

TOTAL SUM OF SQUARES = 0.9867355977049896D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000)00D 04 
F-STATIS{IIC = 0.5821433971075596D 00 
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F-STATISTICS FOR THE 16 COMPONENT 

SUM OF SUMS = 0.6 961387326450431D 02 

SUS OF SUMS**2 = 0.1137157033969463D 03 

SUM OF SQUARE OF SUMS = 0.3597053597878584D 03 

AMCNG GROUP SUM OF SQUARES = 0.3663260306209388D 00 

WITHIN GROUP SUM OF SQUARES = 0.1101186497990677D 03 

TOTAL SUM OF SQUARES = 0.1104849758296886D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.3528623485827816D 00 


F-STATISTICS FOR THE 17 COMPONENT 

SUM GCF SUMS = 0.705285416569i613D 02 

SUM OF SUMS**2 = 0.1548042004914454D 03 

SUM OF SQUARE OF SUMS = 0. 3820324148469269D 03 

AMCNG GROUP SUM OF SQUARES = 0.5041406896350342D 00 

WITHIN GROUP SUM OF SQUARES = 0.1509838763429761D 03 

TOTAL SUM OF SQUARES = 0. 1514880170326111D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D Ou 
F-STATISTIC = 0.3541763825767676D 00 


F-STATISTICS FOR THE 18 COHPONENT 

SUM OF SUMS = 0.7950452893326030D 02 

SUM CF SUMS**2 = 0. 1454127656122291D 03 

SUM OF SQUARE OF SUMS = 0.5157161626920556D 03 

AMNCNG GROUP SUM OF SQUARES = 0.9431815463208067D 00 

WITHIN GROUP SUM OF SQUARES = 0. 1402556039853086D 03 

TOTAL SUM OF SQUARES = 0.14119878553 16294D 03 

DEGREES OF FREEDOM IN NUMEKATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.2148500000000000UD 04 
F-STATiSTIC = 0.7133020797581566D 00 
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F-STATISTICS FOR THE 19 COMPONENT 

SUM OF SUMS = 0.8206200282896035D 02 

SUS OF SUMS*¥*2 = 0. 1642503769158388D 03 

SUM OF SQUARE OF SUMS = 0.5788746875066006D 03 

AMCNG GROUP SUM OF SQUARES = 0.1299298669532475D 01 

WITHIN GROUP SUM OF SQUARES = 0. 1584616300407728D 03 

TOTAL SUM OF SQUARES = 0.1597609287103052D 03 

DEGREES OF FREEDOM IN NOMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STALISTIC = 0.8697276809711277D 00 


F-STATISTICS FOR THE 20 COMPONENT 

SUS OF SUMS = 0.9520048302882186D 02 

SUM OF SUMS**2 = 0.2094086558414699D 03 

SUM OF SQUARE OF SUMS = 0.8713219218815827D 03 

AMCNG GROUP SUM OF SQUARES = 0.2671131239535162D 01 
WITHIN GROUP SUM OF SQUARES = 0.2006954366226541D 03 
TOTAL SUM OF SQUARES = 0. 203366567 8621893D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000009D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000000000uUD 04 
F-STATISTIC = 0.1411744637781583D 01 

P-STATISTICS FOK THE 21 COMPONENT 


“SUM OF SUMS = 0.8994139952150146D 02 


SUS OF SUMS**2 = 0.2335862791765991D 03 

SUM OF SQUARE OF SUMS = 0.7245755635510632D 03 

AMCNG GROUP SUM OF SQUARES = 0.1852785403586404D 01 

WITHIN GROUP SUM OF SQUARES = 0.2263405235410884) 03 

TOTAL SUM OF SQUARES = 0.2281933089446748D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 148500000CQ00000L 904 
F-STATISTIC = 0. 8682828488637504D 00 
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F-STATISTICS FOR THE 22 COMPONENT 

SUS OF SUMS = 0.1035648300735959D 03 

SOM OF SUMS**2 = 0.2414810868436007D 03 

SUS OF SQUARE OF SUMS = 0.9427341493765320D 03 

AMCNG GROUP SUM OF SQUARES = 0.2276892141650119D 01 

WITHIN GROUP SUM OF SQUARES = 0.2320537453498354D 03 

TOTAL SUM OF SQUARES = 0.2343306374914855D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FR “DOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1040764077320930D 01 


F-STATISTICS FOR THE 23 COMPONENT 

SUM OF SUMS = 0.9031234717154746D 02 

SUM OF SUMS**2 = 0.2047768116906234D 03 

SUM OF SQUARE OF SUMS = 0.7149379755971323D 03 

aNCNG GROUP SUM OF SQUARES = 0.171183305488191z2D 01 

WITHIN GROUP SUM OF SQUARES = 0.1976274319346521D 03 

TOTAL SUM OF SQUARES = 0.1993392649895340D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1406000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 
F-STATISTIC = 0.91878 22552243536D 00 


F-STATISTICS FOR THE 24 COMPONENT 

SUB CF SUMS = 0.9171706807335522D 02 

SUM OF SUMS**2 = 0.1900438745495707D 03 

SUS CF SQUARE OF SUMS = 0.8 139562480340389D 03 

AMCNG GROUP SUM OF SQUARES = 0. 2531548763025405D 01 

WITHIN GROUP SUM OF SQUARES = 0.1819043120692303D 03 

TOTAL SUM OF SQUARES = 0. 1844358008322557D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D v2 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1748500000000000uUD 04 
F-STATISTIC = 0. 1476188171339999D 01 
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one F-STATISTICS FOR THE 25 COMPONENT 
SUS OF SUMS = 0.8 440126623444366D 02 


° SUM OF SUMS**2 = 0. 1782254895437147D 03 
a SUS SF SQUARE OF SUMS = 0.6892901073945840D 03 
aH AMCNG GROUP SUM OF SQUARES = 0.2143851912627546D 01 
aN, WITHIN GROUP SUM OF SQUARES = 0.1713325884697689D 03 
ny TOTAL SUM OF SQUARES = 0. 1734764403823964D 03 
x DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
we DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
on F~STATISTIC = 0.1327251499840134D 01 


F~STATISTICS FOR THE 26 COMPONENT 


me SUM OF SUMS = 0.6923579472492147D 02 
ve SUM OF SUMS**2 = 0.1481985058366125D 03 
oe SUM OF SQUARE OF SUMS = 0.4452980334 185913D 03 
, AMCNG GROUP SUM OF SQUARES = 0.1257250153391605D 01 


: WITHIN GROUP SUM OF SQUARES = 0. 1437455255024266D 03 

se TOTAL SUM OF SQUARES = 0. 1450027756558182D 03 

= DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGKEES OF FREEDOM IN DENOMINATOR = 0. 14%85000000000000D 04 

zi F-STATISTIC = 0.9277389287476910D 00 


F~STATISTICS FOR THE 27 COMPONENT 


ie SUM OF SUMS = 0.6366092822129744D 02 

= SUM OF SUMS##2 = 0.1177906819912730D 03 

b SUM OF SQUARE OF SUMS = 0.3935159041535181D 03 

PR AMCNG GROUP SUM OF SQUARES = 0.1233349853537058D 01 
+ WITHIN GROUP SUM OF SQUARES = 0. 1138555229497378D 03 


Oa TOTAL SUM OF SQUARES = 0.1150888728032749D 03 
ror DEGREES OF FREEDOM IN NUMERATOR = 0.172400000000000000D 02 


a DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
a F-STATISIIC = 0.1149027974258136D 01 
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F-STATISTICS FOR THE 28 COMPONENT 

SUM OF SUMS = 0.5387592494868525D 02 

SUM OF SUMS**2 = 0. 1186133442790129D 03 

SUM OF SQUARE OF SUMS = 0.2792595957915495D 03 

AMCNG GROUP SUM OF SQUARES = 0.8575190985312517D 00 

WITHIN GROUP SUM OF SQUARES = 0.1158207483210974D 03 

TOTAL SUM OF SQUARES = 0.1166782674196286D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.14850000000000090D 04 
F-STATISTIC = 0.7853366268738306D 00 


F-STATISTICS FOR THE 29 COMPONENT 

SUS OF SUMS = 0.4536144151906700D 02 

SUM OF SUMS**2 = 0.6564305139058856D 02 

SUM OF SQUARE OF SUMS = 0.1919053830677748D 03 

AMONG GROUP SUM OF SQUARES = 0.5472802462192574D 00 

WITHIN GROUP SUM OF SQUARES = 0.6372399755991081D 02 

TOTAL SUM OF SQUARES = 0.6427127780613007D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 
F-STATISTIC = 0.9 109723144851782D 00 


F-STATISSICS FOR THE 30 COMPONENT 

SUM OF SUMS = 0.3732568029970994D 02 

SUM OF SUMS**2 = 0.5705765228634419D 02 

SUM OF SQUARE OF SUMS = 0.1284820507967475D 03 

AMCNG GROUP SUM OF SQUARES = 0.3560162347433721D 00 

WITHIN GKOUP SUM OF SQUARES = 0.5577283177837672D 02 

TOTAL SUM OF SQUARES = 0.5612884801312009D 02 

DEGAEES OF FREEDOM IN NUMERATOR = 0.14000000060000000 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.6770886363437515D 00 
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F-STATISTICS FOR THE 31 COMPONENT 

SUM OF SUMS = 0.2973876286269883D 02 

SUM OF SUMS**2 = 0.4217833610659068D 02 

SUM OF SQUARE OF SUMS = 0.7726205204918280D 02 

AMONG GROUP SUM OF SQUARES = 0. 18302 45094226048D 00 

@ITHIN GROUP SUM OF SQUARES = 0.4140571558609885D 02 t 
TOTAL SUM OF SQUARES = 0.4158874009552145D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.4688645251806333D 00 


F-STATISTICS FOR THE 32 COMPONENT 

SUM OF SUMS = 0.2142367689889270D 02 

SUM OF SUMS**2 = 0.1414820425416038D 02 

SU“ OF SQUARE OF SUMS = 0.3565723402182890D 02 

AMONG GROUP SUM OF SQUARES = 0.5058971897285643D-01 

WITHIN GROUP SUM OF SQUARES = 0.1379 163191394209D 02 

TCIAL SUM OF SQUARES = 0. 1384222163291455D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 . 
F-STATISTIC = 0.3890854828465455D 00 


F-STATISTICS FOR THE 33 COMPONENT 

SUM OF SUMS = 0.1732004830070321D 02 

SUM OF SUMS*#*2 = 0. 1115582097066481D 02 

SUM OF SQUARE OF SUMS = 0.2312899538600082D 02 

AMCNG GROUP SUM OF SQUARES = 0.3130057176754676D-01 

A4ITHIN GROUP SUM OF SQUARES = 0. 1092453101680480D 02 

TOTAL SUM OF SQUARES = 0. 1095553158857235D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.21485000000000000D 04 
F-STATISTIC = 0.3039120267386334D 00 
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F-STATISTICS OF MEAN VARIANCE, SKEWNESS AND KURTOSIS 


F-STATISTICS FOR THE 1 COMPONENT 

SUB CF SUMS = 0.1111617439422250D 05 

SUM OF SUMS**2 = 0.4319425777749392D 07 

SUS OF SQUARE OF SUMS = 0.1861260733960752D 08 

AMONG GROUP SUM OF SQUARES = 0.3143833877400904D 06 

WITHIN GROUP SUM OF SQUARES = 0.3755407373518861D 07 

TOTAL SUM OF SQUARES = 0.4069790761258952D 07 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.48000000000000000 035 
F~STATISTIC = 0.2870223636519026D 01 


F-STATISTICS FOR THE 2 COMPONENT 

SUM OF SUMS =-0. 3062348822555926D 07 

SUM CF SUMS**2 = 0.7751235739009674D 12 

SUM OF SQUARE OF SUMS = 0.2376364181020636D 13 

AMCNG GROUP SUM OF SQUARES = 0.5306562101878761D 11 

WITHIE GROUP SUM OF SQUARES = 0.7031125381184634D 12 

TOTAL SUM OF SQUARES = 0.7561781591372509D 12 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES CF FREEDOM IN DENOMINATOR = 0.4800000000000000D 03 
F-STATISTIC = 0.2587620620217669D 01 


&-STATISTICS FOR THE 3 COMPONENT 

SUM CF SUMS = 0.1087692759375071D 08 

SUM OF SUMS**2 = 0.1368235571001090D 14 

SUM CF SYUARE OF SUMS = 0.3521258942662316D 14 

AMCNG GROUP SUM OF SQUARES = 0 8280430048704694D 12 

WITHIN GtOUP SUM OF SQUARES = 0.1261530754556778D 14 

ZOTAL SUM OF SQUAKES = 0. 1344335055043825D 14 

DEGREES GF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGKEES OF FREEDOM IN DENOMINATOR = 0.4800400000000000D 03 
F-STATISTIC = 0.22504442145S6077D 01 
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